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ABSTRACT 

We extend our effort to calibrate stellar isochrones in the Johnson-Cousins ( BVIc ) and the 2MASS (JHK S ) 
filter systems based on observations of well-studied open clusters. Using cool main-sequence (MS) stars in 
Praesepe, we define empirical corrections to the lLeieune et al.1 color-effective temperature (T^f) relations down 
to r e ff ~ 3600 K, complementing our previous work based on the Hyades and the Pleiades. We apply em¬ 
pirically corrected isochrones to existing optical and near-infrared photometry of cool (T c ff < 5500 K) and 
metal-rich ([Fe/H]= +0.37) MS stars in NGC 6791. The current methodology relies on an assumption that 
color-T e ff corrections are independent of metallicity, but we find that estimates of color-excess and distance 
from color-magnitude diagrams with different color indices converge on each other at the precisely known 
metallicity of the cluster. Along with a satisfactory agreement with eclipsing binary data in the cluster, we view 
the improved internal consistency as a validation of our calibrated isochrones at super-solar metallicities. For 
very cool stars (T e g- < 4800 K), however, we find that B — V colors of our models are systematically redder than 
the cluster photometry by ~ 0.02 mag. We use color-T^f transformations from the infrared flux method (IRFM) 
and alternative photometry to examine a potential color-scale error in the input cluster photometry. After ex¬ 
cluding B - V photometry of these cool MS stars, we derive E(B — V) = 0.105 ±0.014, [M/H]= +0.42 ± 0.07, 

(m - M )o = 13.04 ± 0.08, and the age of 9.5 ± 0.3 Gyr for NGC 6791. 

Subject headings: Hertzsprung-Russell diagram — open clusters and associations: individual (NGC 6791, 
Hyades, Pleiades, Praesepe) — stars: distances — stars: abundances — stars: evolution 


1. INTRODUCTION 

Although the theory of stellar structure and evolution is 
considered one of the most successful developments in astro¬ 
physics, there still remains a significant mismatch between 
theoretical stellar models and the observed main-sequence 
(MS) of the best studied nearby open clusters. Because MS- 
fitting on color-magnitude diagrams (CMDs) serves as the 
principal method of constructing an accurate distance scale in 
the local universe, identifying and removing underlying sys¬ 
tematic errors in theoretical models not only helps us under¬ 
stand properties of individual stars, but also has a profound 
impact on other fields in astronomy. Since the MS fitting can 
provide information on both distance and foreground redden¬ 
ing, it is an important and useful way of achieving accurate 
luminosity calibration, which is demanded in many astrophys¬ 
ical problems. 

In this series of papers, we have provided an overall as¬ 
sessment of theoretical stellar models and improved the ac¬ 
curacy of models in the Johnson-Cousins and the 2MASS 
(ISkrutskie et alJl200 6) filter systems. We constructed stellar 
modelsjising the Yale Rotating Evolutionary Code (YREC; 
[Sills et aLI l2000t> suite of programs, which provides consis¬ 
tent results with helioseismic observations and solar neu- 


trino fluxes from the Sun (Basue 

al. 2000; Bahcall et al. 

2001; Bahcall & Pinsonneault 2004 

). In Pinsonneault et al. 


(120031 hereafter Paper I), we verified that YREC models are 
in agreement with the masses and luminosities of the well- 
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studied Hyades eclipsing binary system, V818 Tau (= vB 22; 
iTorres & Ribasl2002h . We further demonstrated that our mod¬ 
els provide a good match to spectroscopically determined 
temperatures of individual Hyades members (iPaulson et all 
120031) with good trig onometric parallaxes from Hipparcos 
(Ide Bruiine et all200ll) . 

While YREC models are in satisfactory agreement with the 
best available data on a theoretical plane (i.e., mass, lumi¬ 
nosity, temperature), the match to data becomes poor when 
models are co mpared on an observationa l plane with broad¬ 
band colors. In lPinsonneault et~akl (120041 hereafter Paper II), 
we employed the color-eff ecti ve temperature ( T^ff) relation s of 
lAlonso et al J (1 1 995L fl 996) and lLeieune et al.l (fl997lll998h . to 
transform a theoretically predicted T e ff into observed broad¬ 
band colors. The former color-T'eff relation is based on the 
infrared flux method (IRFM), while the latter is from theoret¬ 
ical spectra with empirical corrections. We found that all of 
these widely used relations fail to reproduce the observed MS 
of the Hyades, with differences in colors as large as 0.1 mag in 
the Johnson-Cousins and 2MASS filter systems. Since YREC 
correctly predicts the mass-luminosity-7’ e ff relation for these 
stars, a natural explanation for the observed offset is a large 
systematic error in the adopted color- T e n relations, at least in 
the parameter space covered by the Hyades members. 

To ease the tension between models and observations, we 
int roduced empirical correctio ns for the color- T e g relations 
of iLeieune et al.1 (119971 [1998) to match photometry of the 
Hyades’ MS (Paper II). This is a much simpler, but more prac¬ 
tical way of overcoming the difficulty than directly examining 
stellar atmosphere models that have large theoretical com¬ 
plexities and uncertainties. Our correction scheme has im¬ 
portant limitations, but we demonstrated that models with the 
Hyades-based color corrections successfully reproduce ob¬ 
served MSs of other nearby, well-studied open clusters (the 
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Pleiades, Praesepe, M67, and NGC 2516) and provide accu¬ 
rate estimates on cluster’s distances^ reddenings, and metal 
abundances from photometry (lAn et al.ll2007bl hereafter Pa¬ 
per III). 

We extended the Hyades-based color- T e ff corrections to hot¬ 
ter and brighter stars using the Pleiades (lAn et al.l2007al . here¬ 
after Paper IV). As described in Paper IV we applied the 
calibrated models to estimate distances and reddenings for 
young Galactic open clusters with Cepheid variables, and 
performed the luminosity calibration of Leavitt’s Cepheid 
period-luminosity (PL) relations. The PL relations con¬ 
structed in this way provided a distance to the active galaxy 
NGC 4258 consistent with itsjnaser-based distance estimate 
(see also lHumphrevs et al.ll2013l . for a recent revision), essen¬ 
tially closing a loop of the distance scale in the local universe. 

Table [H summarizes a library of stellar isochrones with the 
previous and current empirical corrections in this series and 
the accompanying series of papers in the Sloan Digital Sky 
Survey tSDSS: lAlrneTaiml ugriz system. The Hyades- 
based corrections in Paper II are valid only at My < 8 or 
7)if > 4000 K because of the magnitude limit of Hipparcos 
parallaxes for faint Hyades stars dPerryman et al.llT~998i) . The 
Pleiades could not be used to extend the calibration because 
its lower mass stars have not yet arrived on the MS, and K 
dwarfs in the cluster seem anomalously blue in broadband 
colors (IStauffer et alJl2003l see also Paper III). In §|2]we de¬ 
scribe an extension of empirical color calibration to cooler 
stars (Teff > 3600 K) using low-mass stars in Praesepe. This 
complements our earlier color corrections to YREC, cover¬ 
ing a wider range of effective temperature of stars (3600 K 
< T efs < 12000 K). 

In § [3] through § 0 we describe the application of the 
calibrated set of isochrones to NGC 6791, a benchmark 
metal-rich system in the Galaxy, taking advantage of high- 
quality optical and infra red (IR) photometry (Stetso rTet al.l 
120031 ICarnev et al.l2005l) and precise measureme nts of eclips¬ 
ing binary systems (IBrogaard et al.l 1201 ll 120121) . There are 
a considerable number of studies on astrophysical parame¬ 
ters of the cluster in the literature. Early photometric stud¬ 
ies suggested that_the cluster’s metallicity is_near solar (e.g., 
lAnthonv-Twarog & Twarogll 1 9851 [Geisler et aklll99lL among 
others), while more recent studies based on low-resolution 
spectra ind icated [Fe/H]~+0.2 (e.g.. iGarnavich et al.lll994l 
iFriel et al.l [20021 among others). However, the most re¬ 
cent efforts using high-resolution spectroscopy consistently 
found [Fe/H ]~ +0.4. Table 0 summarizes high-resolution 
spectroscopic metallicity estimates and their errors in the 
literature. Most of the studies found an ~ 0.1 dex er¬ 
ror in the mean [Fe/H] of the cluster, with additional sys- 
tematic errors of ~ 0.1 dex as shown in the last column 
dGratton et al.ll20 06: Orig lia et~aDl2006llBrogaard et al.ll20 1 ll 

Boesg aard et al.ll2015l) . Throughout the paper we refer to an 
unweighted mean and a standard deviation ([Fe/H]= +0.37 ± 
0.07) of these measurements as a spectroscopic metallicity of 
the cluster. 

The o-element abundances of star s in NGC 6791 are 
essentially solar JC arraro et al.l 120061 lOriglia et ali 120061: 
Carretta et al.l 120071 iBrogaard et al.l 1201 It Boesgaa rd et al.l 
20151) . Althou gh there is some evidence that [C/Fe] i s below 
the solar value (IGratton et al .120061 lOriglia et al.120 06). a sub- 
solar fO/Fel among cluster members is still in dispute (e.g. , 
Gratton et al.l l2006t lOriglia et al.l 120061 iGeisler et al.l 120121 

iBoesgaard et al.ll2015llCunha et al.l2015l) . In any case, the ef¬ 

fect of CNO abundances is likely to be small in most parts 


of MS compared to a change in the bulk metallicity (e.g., 
IBrogaard et al.ll20l2h . We assumed a scaled solar abundance 
of the cluster, and used [M/H] and [Fe/H] interchangeably to 
indicate the bulk metallicity of its stars. 

While there has been a convergence in the cluster’s metal¬ 
licity, the foreground reddening is not well constrained. 
A group of studies found that the mean foreground red¬ 
dening of NGC 6791 is around E(B - V) = 0.10 (e.g., 
IStetson et al] 120031 ICarraro et al] 120061 among others). On 
the other hand, there is a broad distribution of E(B^_V) esti¬ 
mates centered at E(B - V^ 0.15 (e.g jKaluznv & Rucinskil 
119951 lAnthonv-Twarog et al.1120071 amo ng others). Mean¬ 
while, dust emission maps (lSchlegel et al l 119981) set the in¬ 
tegrated line-of-sight extinction toward NGC 6791 (\b\ ss 
11°), which is E(B - V) « 0.1 33 after a 14% dow nward 
revision in the reddening scale (ISchlaflv et al.l 12010 ). Re¬ 
cent estimates of the true distance modulus of the clus¬ 
ter range from (m — M)p ~ 12.9 to (m - M) n ~ 13.1 (e.g., 
Salari s et al.ll2004l ICarnev et al.ll2005c lAnthonv-Twarog et al.l 

20071 IBrogaard et al.ll201 ll among others). A dispersion in 
these distance estimates is partly due to uncertainty in redden¬ 
ing, since A E(B - V) = 0.05 is translated into A (m - M) o « 
0.12 in the MS fitting. When it comes to a determination 
of absolute magnitudes (My) of stars in the cluster, errors 
in My from distance [a(My) = <T( m -M) 0 ~ 0.12] and extinc¬ 
tion errors [er (My) = Rv<je(b-v ) ~ 0.16], both of which are in¬ 
duced by the error in reddening [<7e(b-v) = 0.05], are added 
together (positively correlated) to give a(My) ps 0.28 mag, 
where Ry = 3.26 is the ratio of total to selective extinction. 
Certainly, more work is needed to shrink the size of the error 
in the cluster’s E(B - V). 

The primary goal of this paper is to extend empirical color- 
T e ff corrections using cool MS stars in Praesepe and to validate 
the calibration at super-solar metallicity using multicolor pho¬ 
tometry of NGC 6791. In the following section, we describe 
the extension of color-TAr corrections. In § [3] we describe 
how we collected and combined optical and IR photometry 
of NGC 6791. Before we apply color-calibrated models to 
the photometric data, we compare, in § [4] theoretical stellar 
models against precise measurements of eclipsing binary sys¬ 
tems in NGC 6791, and show that our interior models are in 
satisfactory agreement with the data with reasonable assump¬ 
tions on r e ff scale and helium abundances. In §0we conduct 
tests on the calibrated models using photometry of NGC 6791 
with the cluster’s well-known spectroscopic metallicity. We 
also report our best-fitting photometric estimates of the clus¬ 
ter’s reddening, distance, metallicity, and age, along with a 
thorough evaluation of their errors. These parameters will be 
used in a subsequent paper (D. An et al. 2015, in prepara¬ 
tion) to derive accurate luminosities and masses of red-clump 
stars, providing useful information on the amount of mass loss 
along the red giant branch in NGC 6791. 

2. EXTENSION OF EMPIRICAL COLOR CALIBRATION 

The general scheme of our calibration procedure is de¬ 
scribed in Paper II, where we forced a match between com¬ 
puted and observed cluster sequences on CMDs by adjusting 
color-T e ff relations in Leieune et all (fl997Ll!998h . Empirically 
corrected isochrones constructed in this way were designed to 
match a MS of a calibrating cluster system. In Paper II we 
used the Hyades in the calibration, but the Hyades-based cor¬ 
rections were defined at My < 8 or 7Ar > 4000 K because 
of the lack of faint Hyades members with good Hipparcos 
parallaxes dPerryman et al JM 9981) . Our new empirical color 
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TABLE 1 

A Library of Empirical Calibrated Isochrones 


Reference 

Interior 

Models 

Base Color-r e ff 
Relations 

Filter 

Passbands 

Calibrating 

Clusters 

Valid Range 
in r eff 

Papers II and III 
Paper IV 

This work 

YREC (OPAL) 
YREC (OPAL) 
YREC (OP) 

Leieune et al. 
Leieune et al. 
Leieune et al. 

BVI C K S 

BVI C JHK S 

BVI C JHK S 

Hyades 

Hyades, Pleiades 

Hyades, Pleiades, Praesepe 

4,000-8,200 K 
4,000-12,000 K 
3,600-12,000 K 

An et al. (2009b) 
An et al. (2013) 

YREC (OP) 
YREC (OP) 

MARCS 

MARCS 

ugriz 

ugriz 

~M57 

Primarily M67 and M92 

4,000-6,000 K 
4,000-7,000 K 


TABLE 2 

Metallicity of NGC 6791 from High-Resolution 
Spectroscopy 


Reference 

[Fe/H] 

cr (internal) 

cr (systematic) 

Gratton et al. (2006) 

+0.47 

0.04 

0.08 

Carraro et al. (2006) 

+0.39 

0.01 


Origlia et al. (2006) 

+0.35 

0.02 

<0.1 

Brogaard et al. (2011) 

+0.29 

0.03 

0.07 

Geisler et al. (20 i 2) 

+0.42 

0.01 


Cunha et al. (2015) 

+0.34 

0.06 


Boesgaard et al. (2015) 

+0.30 

0.02 

$0.1 

Unweighted Mean 

+0.37 

(LOT 3 



a Unweighted standard deviation. 


correction procedure described below extends the calibration 
further down below this limit by employing CMDs of Prae- 
sepe. Praesepe has a similar age (550 Myrf] and metallic¬ 
ity ([Fe/H]=+0.14, see below) as the Hyades, but exhibits 
a negligible depth effect at its intermediate distance from 
the Sun [(m — M )o = 6.33 ± 0.04; Paper III]. At the hot end 
(?eff > 8500 K), the Hyades-based calibration was limited 
due to the relatively old age of the cluster, so we used young 
Pleiades stars to extend color corrections to hotter stars (Pa¬ 
per IV). Our revised color calibrations presented in this paper 
are based on all of these benchmark open clusters, and super¬ 
sede the previous color corrections. Luminosities of the mod¬ 
els in these studies are ultimately tied to the Hipparcos^ paral¬ 
laxes to the Hyades (m - M )o = 3.33 ± 0.01 (Pe rryman et all 
U998h . 

2.1. Base Models 

The isochrone method assumes that the luminosity-7]^ re¬ 
lationship in the models is correct, and maps temperature 
onto color via the inferred distance, magnitudes, and bolo- 
metric corrections. In this paper, we used an updated ver¬ 
sion of theoretical isochrones generated using the YREC 
stellar evolutionary code (F. Delahaye & M. H. Pinson- 
neault, in preparation). These theoretical models depend 
on the choice of input physics, in particular the opacity ta¬ 
bles and chemical mixture. In this work we have switched 
to using the atomic op acity tables provided by the Opac¬ 
ity Project (OP: iBadnell etakll2005|l from the OPAL tables 
(llglesias & Rogers! [ 19961) . which slightly alters the underly¬ 
ing isochrones in the theoretical plane (see the second col¬ 
umn in Table |T}. The current models supercede the prior 
ones. Detailed information on the OP-b_ased YREC models 
can be found in iDelahave & Pinsonneauli (2 0051 2006 ): see 

4 The age of the H yades is about 625 Myr based on models with convective 
overshooting IPerrvman et all 1 998 1, but a difference in age has little impact 
for low-mass MS stars. 


also lvan Saders & Pinsonneaulti (120121) for an updated discus¬ 
sion of input physics. Even though the opacity approaches are 
different, they have similar Rosseland means for solar condi¬ 
tions and thus yield similar solar calibrations. In the temper¬ 
ature range where these models were compared to MS stars 
in NGC 6791 (3600 K < r eff < 5600 K), the biggest differ¬ 
ence (ATetf ~ 100 K) is found at 7)n ~ 4600 K in the sense 
that the OP-based models predict hotter temperatures. Below 
we repeated our earlier color calibrations based on the Hyades 
(Paper II) and the Pleiades (Paper IV) using OP-based models 
for an internal consistency. 

The remaining paramet ers in the models, includ ing the 
initial chemical mixture (iGrevesse & Sauvall 119981) . mix¬ 
ing length (a = 1.72), and helium enrichment parameter 
(AT/AZ = 1.2), remained the same as in our previous 
papers of this series. The helium enrichment parameter 
was d erived from the p r imordial helium abund ance (Y p = 
0.245; iBono et aD 120021 : iThuan & Izotovl 120021 and refer¬ 
ences therein) and solar models that lack microscopic diffu¬ 
sion. As shown in Paper I, models without diffusion cor¬ 
rectly predicted luminosities of the vB 22 eclipsing binary 
system in the Hyades. However, the microscopic diffusion 
of heavy elements and helium could be important at the age 
of NGC 6791. Nevertheless, we neglected microscopic dif¬ 
fusion in our analysis, because there is as yet no strong ob¬ 
servational evidence of such effects in the cluster (see discus¬ 
sions in lBrogaard et alJ[2012l) . In any case, models with mi¬ 
croscopic diffusion are only few tens of Kelvin hotter in most 
parts of the cluster’s MS (e.g. JChabover et ahbOOll) . while a 
change of heavy metals in the surface of stars is likely to be 
small (A[Fe/H] < 0.05) along the MS fe.g.. iBrogaard et ail 
12012L see their Figure 7). 

The YREC evolutionary tracks were interpolated to gen¬ 
erate theoretical isochrones at stellar ages from 1 Gyr to 
16 Gyr. Models were constructed over-0.3 < [Fe/H] < +0.5 
in 0.1 dex increments and at [Fe/H] =+0.75. Luminosities 
and T e ff were initially converted to V,B — V,V — Ic,V — K s , 
J - K s , and 77 - K s in the Johnson-C o usins and 2MASS sys¬ 
tems using the iLeieune et al.l (119971119981) color-7^jf table, 
where “merged” color-7[n relations of their original tables 
(see Paper IV) were used in order to produce smooth base 
isochrones. 

In Figures Q] and [2] isochrones are compared to CMDs of 
Praesepe with B - V, V - Ic, V - K s , J - K s , and H - K s as the 
color index (hereafter BV, Vic , VK S , JK S , and HK S CMDs, re¬ 
spectively), where we used the same cluster photometry as in 
Paper III. Here, we present CMDs of Praesepe to highlight the 
extension of color- T e s calibration in the lower MS. Interested 
readers are referred to papers in this series for CMDs of the 
Hyades and the Pleiades used in the following analysis. 

The plus symbols in Figures Q] and [2]represent cluster mem¬ 
bers after excluding previously known binaries. Because of 
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Fig. 1.— CMDs of Praesepe. Plus signs are cluster members after excluding previously known binaries. Open circles are those remained from a statistical 
filtering routine, which was design ed to ide ntify and remov e unresolved binaries and/or blended sources independently of models. Dotted lines are 550 Myr old 
models at [Fe/H]= +0.14 using ILeieune et all d 1997L119981) color-relations, and thick solid lines represent models with empirical corrections derived in this 
paper. For the purpose of comparison, calibrated isochrones at 650 Myr are additionally shown as a thin solid line. Isochrones are placed at (m — M )o = 6.33 
assuming E(B -V) = 0.006. 



J - K H - K 

S S 


Fig. 2.— Same as in FigureQ] but CMDs with J — K s and H — K s indices. 
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extensive membership and binarity data in the literature (see 
references in Paper III), most parts of the cluster’s CMDs are 
relatively free of cluster binaries and foreground/background 
stars except at the bottom of MS. We employed a photometric 
filtering routine that utilized photometric data in various color 
indices to identify and reject likely cluster binaries and/or 
background stars independently of theoretical isochrones (see 
Paper III). Open circles are those remained after the photo¬ 
metric filtering, and were used in the following calibration 
process. 

Th e dotted lines i n Figures 0 and 0 are YREC isochrones 
with iLeieune et all (I1997L fl998 ) colors. The solid lines are 
models with our empirical color corrections that are described 
below in detail. We assumed the same cluster age with no 
core overshoot (550 Myr) as the Hyades (e.g., iMermilliodl 
ll98!l) . and adopted the best estimates on the cluster metallic - 
ity ([Fe/H] = +0.14 ± 0.02) and foreground reddening [E(B - 
V ) = 0.006 ± 0.002] in Paper III. More specifically, the clus¬ 
ter metallicity is the weighted mean of our spectroscopic 
([Fe/H] = +0.11 ±0.03) and photometrically derived metal¬ 
licity ([M/H] = +0.20±0.04). The former is based on Fe I and 
Fe II line measurements of four MS stars in the cluster, using 
the Magellan Inamori Kyocera Echelle (MIKE) spectrograph 
(Bernstein et al.ll2003lf on the Magellan 6.5-m Clay telescope. 
The photometric metallicity was determined using isochrones 
with the Hyades-based color calibration. Our adopted metal¬ 
licity o f Praesepe is consistent w ith more recent determina¬ 
tions inlC arrera & Pancinol(l201 lL [Fe/H] =+0.16±0.05) and 
iBoesgaard et all (I2013L fFe/H] = + 0.12 ± 0.04), but is lower 
than the value in lPaceet all (120081 [Fe/H] = +0.27±0.10). 
Our adopted cluster reddening is the average of E(B — V ) es¬ 
timates in the literature and our own photometric estimates 
using calibrated isochrones in Paper III. We also took the 
cluster’s distance [(m - M )o = 6.33 ± 0.04] from Paper III, as 
derived using the Hyades-based calibrated models with the 
above mean metallicity and reddening. In Figures 0 and [2} 
650 Myr old models with color corrections are shown in thin 
solid lines to illustrate the age dependence of models. 

Throughout the paper, we adopted extinction laws in Pa¬ 
per IV (equations 2-7). Briefly, these e quati ons have color- 
excess ratios for zero-color stars from ICardelli et al.l ( 1989 ). 
i.e., E(V - I C )/E(B -V)= 1.30, E(V - K S )/E(B - V) = 2.88, 
E(J - K S )/E(B -V) = 0.56, and E(H - K S )/E(B - V) = 0.20 
with Ry = 3.26 for the ratio of total-to-selective extinction. 
They include color terms in the color-excess ratios and Ry 
from lBessell et al.l ( 1998 ) to take into account effective wave¬ 
length shifts in broadband filters for cooler stars. The cluster’s 
E(B -V) in the following analysis refers to a value for zero- 
color stars. 

2.2. Empirical Color Calibration 

Figures 00 show differences between theoretical and ob¬ 
served colors of individual stars in the Hyades, the Pleiades, 
and Praesepe, respectively. Color differences (in the sense 
of observed colors minus model values) were computed us¬ 
ing models without color calibration (e.g., dotted lines in Fig¬ 
ures Q] and [2j at constant My in B — V, V — Ic, and V - K s 
colors, and at constant M K in J - K s and H — K s . The tem¬ 
peratures on the abscissa are 7/n inferred from My or Mk in 
the models. We adopted the same cluster parameters (age, 
metallicity, reddening, and distance) of the Hyades and the 
Pleiades as in the previous papers of this series: 550 Myr, 
[Fe/H] = +0.13 ± 0.01, E(B - V) = 0.000 ± 0.002, (m - M) 0 = 
3.33 ±0.01 for the Hyades, and 100 Myr, [Fe/H] =+0.04± 


0.02, E(B - V) = 0.032 ± 0.003, (m - M) 0 = 5.63 ± 0.02 for 
the Pleiades. The Pleiades’ distance is the average geomet¬ 
ric distance to the cluster (see Paper III), and is consistent 
with the most recent determination of trigonometric paral¬ 
lax using very long baseline radio interferometry dMelis et al.l 
12014 ). Single cluster members of the Hyades after exclud¬ 
ing known binaries are shown in Figure 0 Photometrically 
selected single stars are shown for the Pleiades in Figure [4] 
after rejecting known cluster binaries and non-members (see 
references in Paper III). Cool MS stars in the Pleiades with 
T e ff 5000 K are not displayed, because the young age of the 
cluster makes broadband colors of its low- mass s tars deviate 
from those observed in older systems (Stauffer et al.l2003l see 
also Paper III). 

As shown in Figures 00 color deviations of the models 
from the three cluster systems look similar, suggesting a sys¬ 
tematic orig in of these offsets. We defi ned empirical correc¬ 
tions for the ILeieune et al.l ( 1997 , 1998 ) color-7]^ relations as 
these systematic color residuals, and are shown as solid lines. 
The detailed steps in constructing these curves are described 
below. 

The vertical dotted lines in Figures 00 show 7/n ranges, 
in which each of the calibrating cluster systems was em¬ 
ployed in the empirical color-7±f corrections. At 4500 K 

< 7 d -f < 7500 K, we employed the Hyades, and computed 
moving averages of color differences using 4-8 data points 
in each moving window, with a 50% overlap between adja¬ 
cent windows. The upper 7/n limit was set to 7/n = 5500 K 
in JK S and HK S CMDs, because of large photometric errors 
of saturated stars in 2MASS. We used 7-point moving aver¬ 
ages of the Pleiades stars with 50% overlaps to define color 
corrections at 7/n > 8500 K in BV, Vic , and VK S CMDs, 
and at 7/n > 6800 K in the JK S and HK S CMDs. At 7500 K 

< 7)if < 8500 K, the Pleiades stars show smaller color devia¬ 
tions than those in the Hyades, which may be due to a younger 
age and/or a lower metallicity of the Pleiades than the Hyades, 
but other error sources, such as a photometric zero-point er¬ 
ror, may also be responsible for the discrepancy. We used a 
linear ramp to bridge over the gap (7500 K < T e ff < 8500 K 
or 5500 K < 7 dl < 6800 K), where the ramped curves simply 
represent intermediate color deviations from the two clusters. 

The systematic trend observed in Praesepe (Figure [5]) is 
in good agreement with that of the Hyades (Figure 0 over 
4000 K < 7/ff < 8000 K. Given the similar age and metallic¬ 
ity of the two systems, the observed similarities strongly sug¬ 
gest that the color residuals are systematic in nature. We em¬ 
ployed Praesepe to define color corrections at 7/n < 4200 K 
with 4-point moving averages and 50% overlaps between ad¬ 
jacent moving windows. A linear ramp was used over 4200 K 

< 7jff < 4500 K to combine color corrections from Praesepe 
and the Hyades. 

We smoothed the moving-averaged points by taking a lin¬ 
ear interpolation of five neighboring points. We used three 
moving-averaged points in the smoothing in the upper and 
the lower ends of T e g. Our final set of empirical corrections 
are shown as solid lines in Figures 00 and are tabulated in 
Table 0 We assumed that our empirical corrections are in¬ 
dependent of metallicity and age (see below), and applied to 
all isochrone s generated using YREC and the Leieune et al. 
(119971 fl998h color transformations Q 

In Paper III we already demonstrated that isochrones with 
the Hyades-based corrections provide excellent matches to the 

5 Available at http://home.ewha.ac.krrdeokkeun/astro/isochrone.html 
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Fig. 3.— Color differences between models using ILeieune et all dl997LlT998h relations and single cluster members in the Hyades with accurate parallaxes. A 
550 Myr old model with [Fe/H]= +0.13 was used at (m — M )o = 3.33 without foreground dust extinction correction. Color differences were computed at constant 
My for B -V,V — Ic, and V — K s , and those at constant for J — K s and H — K s , in the sense of observed data minus model values. Solid lines are empirical 
color-r e ff corrections, or moving-averaged points of the systematic color difference, constructed from the three calibrating cluster systems: the Hyades (4500 K 
< 7rff < 7500 K, or T e ff < 5500 K in / — Ks and H — K s ), the Pleiades (F e ff > 8500 K or 6800 K in / — Ks and H — Ks; see Figure|4j, and Praesepe (7 e ff < 4200 K; 
see Figure [5)• The vertical dashed lines indicate these 7 e ff divisions. 


TABLE 3 

Empirical Color Corrections 


Rtf 

(K) 

A (B - V) 
(mag) 

MV - I c ) 

(mag) 

Mv - Ks) 
(mag) 

A (/-Ks) 
(mag) 

A (H - K s ) 
(mag) 

3600 

0.0877 

0.2857 

0.3947 

0.0810 

0.0200 

3650 

0.0877 

0.2834 

0.3947 

0.0810 

0.0200 

3700 

0.0877 

0.2546 

0.3817 

0.0845 

0.0241 

3750 

0.0868 

0.2231 

0.3548 

0.0905 

0.0291 

3800 

0.0856 

0.1948 

0.3256 

0.0955 

0.0311 


NOTE. — Only a portion of this table is shown here to demonstrate its form and 
content. A machine- readable version of the full table is available. 


MS of nearby open clusters, including the Pleiades (see Fig¬ 
ure [4]), Praesepe (see Figure 0, M67 ([Fe/H]=+0.00±0.01, 
4 Gyr), and NGC 2516 ([Fe/H]= -0.04 ±0.05, 140 Myr). 
This implies that our empirical color corrections are primar¬ 
ily a function of T e ff in the metallicity range covered by these 
systems (—0.04 < [Fe/H] < +0.14). However, it is necessary 
to check and validate models beyond the above metallicity 
range. The main goal of this work is to perform a stringent 
test of empirical corrections at the high metallicity end using 


NGC 6791. 

In Figures Q] and O isochrones incorporating these correc¬ 
tions are shown as solid lines on the Praesepe CMDs. Prae¬ 
sepe was solely used to define color offsets at r e ff < 4200 K, 
which correspond to B — V « 1.3, V — Ic ~ 1.6, V — K s « 3.3, 
and / - K s ~ 0.8. Colors are based only on the Hyades stars 
at 0.3 <B-V< 1.2, 0.4 < V - I c < 1.3, 0.8 < V-K s <2.8, 
and 0.4 < J - K s < 0.7. Model colors near the cluster’s MS 
turn-off hinge on calibrations from both the Hyades and the 
Pleiades. As noted above, the two systems predict slightly dif¬ 
ferent color corrections in this region, which partly explains 
why our 550 Myr old model shows a color offset near the clus¬ 
ter’s turn-off. However, this offset should not affect our results 
in the following analysis, because old systems like NGC 6791 
have MS stars cooler than r e ff ~ 7000 K. 

2.3. Comparison with 1RFM Color-T e ff Relations 

The IRFM has been widely used to obtain temperatures 
of stars based on photometric observations, with little de¬ 
pendence on theoretical stellar at mosphere models (see ref¬ 
erences in lCasagrande et akll20l O). Considering that 7)n is a 
defined quantity from L = 4~R 2 irT^ t> the IRFM temperatures 
































































Distances from Main-Sequence Fitting. V 


7 



T eff (K) 



4.1 4.0 3.9 3.8 3.7 3.6 4.1 4.0 3.9 3.8 3.7 3.6 

log T eff (from M v ) log T eff (from M K ) 


Fig. 4.— Same as in Figure [5] but showing color differences between model and observed colors of single cluster members in the Pleiades. Stars below 
Teff ~ 5000 K are not displayed because of anomalous colors related to the young age of the cluster. A 100 Myr old model at [Fe/H]= +0.04 was used with 
(m - M) 0 = 5.63 and E{B - V) = 0.032. 


are arguably closer to a fundamental scale than spectroscopic 
temperatures, either based on an excitation equilibrium of Fe I 
lines or Balmer line profiles. On the other hand, a fundamen¬ 
tal scale can be obtained independently using observations of 
well-studied cluster systems with accurate distance, redden¬ 
ing, and metallicity, as described in this series of papers. 

The two independent approaches yield the same colors of 
the Sun within errors. At the age of the Sun (4.57 Gyr), 
colors from our calibrated solar metallicity isochrone are 
B -V = 0.651, V-I c = 0.704, V-K s = 1.563, and J-K s = 
0.336 at r R ff.r?i = 5777 K. The IRFM colors of the Sun in 
ICasagrande et al.1 (1201 Ol) are B - V = 0.641 ± 0.024 ± 0.004, 
V -I c = 0.690 ±0.016 ±0.004, V - K s = 1.544 ± 0.018 ± 
0.010, and J - K s = 0.362 ± 0.029 ± 0.003, where errors in 
each quantity represent a random and a systematic error, re¬ 
spectively. 

In the left panels of Figure[6]our calibrated color-T -it rela- 
tions a re compared to the IRFM relations in Casagran de et al.1 
(2010) over a wide range of T e ff in B - V, V — Ic, V - K„, 
and J - K s , respectively. In these comparisons, IRFM col¬ 
ors were estimated directly from T e g in the YREC isochrones. 
We used 9 Gyr old models, but the comparisons in Fig¬ 
ures [6] largely represent differences between the two color- 
T e ff relations, independently of the adopted age of models. 
Color differences are shown at [Fe/H]= -0.3 (blue solid line). 


[Fe/H]= 0.0 (black dashed line), and [Fe/H]= +0.4 (red solid 
line). Since the IRFM color-T e ff relation in V — Ic is de¬ 
fined up to [Fe/H]= +0.3, the V — Ic comparison is shown at 
[Fe/H]= +0.3 by the red solid line. 

As seen in these comparisons, IRFM relations show better 
agreement with our empirically corrected c olor-7jn relations 
than the original lLeieune et al.1 (119971 fl9 98) relations. In par¬ 
ticular, color differences are significantly reduced to ~ 0.02- 
0.04 mag for cool MS stars (T C ff < 5000 K). Both IRFM and 
the cluster-based approaches are limited in the super metal- 
rich regime by a small number of calibration stars, but the 
agreement in Figure [6] shows that the two independent color- 
7^ff relations are now converging on a consistent temperature 
scale. 

The T e ff differences between the IRFM and the calibrated 
stellar isochrones are shown in Figure [7] in each color index. 
The temperatures from these two approaches are nearly on the 
same scale, although isochrones in / — K s show systematically 
higher temperatures at all colors. Nevertheless, there still re¬ 
main small-scale structures in the color difference between 
IRFM and the calibrated isochrones, which may indicate er¬ 
rors in either of these two approaches. Below we directly 
compare isochrones with an observed MS of NGC 6791 to 
test the accuracy of color-T e ff relations (§0. 

In Figures [6] and [7] the sets of lines with different metal- 
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Fig. 5.— Same as in Figure[3] but showing color differences between model and observed colors of single cluster members in Praesepe. A 550 Myr old model 
at [Fe/H]= +0.14 was used with (m - M )o = 6.33 and E(B - V) = 0.006. 

licities are not identical, indicating that our model and IRFM 
relations have different degrees of sensitivity to metallicity. 

Since our empirical corrections do not have an explicit de¬ 
pendence on metallicity, the differences show that the metal¬ 
licity dependence of the iLeieune et al.1 (119971 179981) relations 
differs from that of IRFM colors. The differences are rela¬ 
tively small, but they are of the order of systematic color or 
T e ff differences for a given metallicity. Photometric estimates 
of cluster parameters (distance, reddening, and metallicity) 
depend on the adopted sensitivity to metallicity of models, 
so these two systems should give slightly different photomet¬ 
ric solutions. In §[5]we discuss the metallicity dependence of 
models using NGC 6791. 

Since IRFM relations in lCasagrande et ahl (120101) treat both 
MS and subgiant stars equally without involving gravity 
terms, we did not include color or T e a comparisons for sub¬ 
giant stars in Figures [6] and [7] The effect of surface gravity 
on stellar colors is generally milder than that of a t empe rature 
or a metallicity change. In ILeieune et alJ d 19971 Il998t) . this 
amounts to ~ 0.01-0.02 mag between MS and subgiant in 
the above broadband colors. In this paper, we derived empir¬ 
ical corrections based on MS stars, and applied them to sub¬ 
giant branches in the models with lLejeune et al.! ( 1997 . 1998 5 
colors. Empirical correction terms that explicitly depend on 
gravity are expected to be less than ~ 0.01 mag, if any. Since 


the following analysis is mainly concerned with MS stars in 
NGC 6791, systematic errors from surface gravities can be 
neglected. 

3. PHOTOMETRY OF NGC 6791 

In our analysis, we utilized optical and near-IR photometry 
of NGC 6791 in the literature. Figure [8] shows the cluster’s 
CMDs, where BVIc photometry originally from lStetson et al. 
(2003) is combined w ith /('-band magnitudes in ICarnev et al. 
(2005). IStetsonl (120051) found color-scale errors in his standard 
star database, and su bsequently repo rted the update of the 
cluster photometry in iBrogaard et all ( 2012 ) with the newer 
calibration. We took the updated photometry table from 
Stetson’s homogeneous photometry websittQ and utilized it 
throughout the anal ysis. In this paper, we keep the orig¬ 
inal reference (IStetson et al.l 120031) to indicate the updated 
photometric table. Taking the same selection criteria as in 
IStetson et al. ( 200 3) for good photometry, we applied cuts 
based on sharp, separation indices, and photomet¬ 
ric errors. We assigned minimum errors of 0.002 mag in V, 
B -V, and V - lc ■ We combined the optical photometry with 
near IR photometry after transforming J- and /C-band mea¬ 
surements bv iCarnev et al. ( 2005 ) in the CIT system into JK S 
in 2MASS, using a color transformation equation found in 

6 http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/STETSON/ 
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6000 5500 5000 4500 6000 5500 5000 4500 

Jeff (YREC) T ef f (YREC) 

Fig. 6.— Comparisons between IRFM JCasagrande et alJI201QD colors and those from 9 Gyr old m odels in B - V, V — Ic, V — K s , and J — K s , in the sense 
of isochrone colors minus IRFM values. Color differences from calibrated models and those using ILeieune et ahl 119971 fT998i) colors are displayed in the left 
and right panels, respectively. Comparisons are shown at three different metallicities: [Fe/H]= —0.3 (blue solid line), [Fe/H]= 0.0 (black dashed line), and 
[Fe/H]= +0.4 (red solid line; [Fe/H]= +0.3 in AV - 7 C ). 


I Carpenter! (1200 ll) FI The sloping dotted line in t he VK S CMP 
in Figure [8] represents the completeness limit in iCarnev et al.l 
(120051) at K = 16.5 mag, which is ~ 2 mag deeper than the 
99.9% completeness limit in 2MASS (K s ~ 14.3 mag). Al¬ 
though the VK S CMD is limited to the upper part of the clus¬ 
ter’s MS, the observed sequence provides usefu l constraints 
on the cluster’s reddening and metallicity (§ 15.11 ). 

Unresolved cluster binaries (or photometric blends) and 
cluster non-members are sources of bias in the MS fitting. 
There exists useful information on the cluster membership 
from proper-motion measurements dPlatais et al.l 1201 ll) . but 
there is no detailed census taken on binaries along the cluster’s 
MS. In order to identify and remove outliers from the cluster 
MS, we applied a photometric filtering routine (Paper III) to 
the BV , Vic, and VK S CMDs. This step was done indepen¬ 
dently of the models. Stars tagged as outliers in one of the 
CMDs were rejected from all of the CMDs in the above fil¬ 
tering procedure. Black cross points in Figure |8]are those re- 

7 Updated color transformations can be found at 
http://www.astro.caltech.edu/~jmc/2mass/v3/transformations/ 


maining after the photometric filtering, with a final \ 2 thresh¬ 
old value corresponding to 2.5a (see Paper III). As seen in 
Figure [8] the photometric filtering routine provides a robust 
detection and removal of binaries and background stars, even 
without cluster membership information. On the other hand, 
lower mass-ratio binaries could still remain after the filtering 
process, because they lie near a single-star MS. We included 
a systematic error induced by these leftovers in th e tot al er¬ 
ror budget of the cluster parameter estimates (see § 15.61) . Be¬ 
cause the filtering routine was designed to select single stars 
on the MS, we rejected outliers by hand near MS turn-off 
(V < 17.6 mag) and on the subgiant branch. 

Open circles in Figure [8] represent single-star cluster se¬ 
quences that were derived using the above remaining stars. At 
V > 17.5 mag, we computed the median color of these stars 
in each magnitude bin having AV = 0.1 mag. We constructed 
a cluster sequence near MS turn-off by picking a number den¬ 
sity peak along the observed sequence. In the subgiant region, 
we computed a median V magnitude in each color bin with a 
bin size of 0.03, 0.03, and 0.08 mag in the BV, Vic, and VK S 
CMDs, respectively. These cluster sequences were used in a 
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6000 5500 5000 4500 
T e „ (YREC) 

Fig. 7.— Same as in Figure [6] but showing T e ff differences between cali¬ 
brated isochrones and IRFM relations, in the sense of the former minus latter. 
The dotted lines indicate AT e ff = 0, ±100 K. 

joint determinatio n of the cluster’s distance, metallicity, red¬ 
dening, and age ($ 15.5k 

4. TEST OF MODELS WITH ECLIPSING BINARIES IN 
NGC 6791 

In Paper I we compared the YREC models with eclipsing 
binary data for the well-studied system (vB 22) in the Hyades, 
from which we found a reasonably good agreement of theoret¬ 
ical models with observed mass-luminosity and mass-radius 
relations. This result provided a basis for the empirical correc¬ 
tions to Teff-to-color transformations, which improved a match 
between theoretical and observed colors of stars in nearby 
clusters. However, NGC 6791 is a more metal-rich system 
than the Hyades and the other nearby clusters used in our cal¬ 
ibration. Since systematic errors in the interior models can be 
propagated into errors in the stellar colors and magnitudes, a 
verification of our empirical corrections should be preceded 
by a test of i nterior models a t t he me t allicity of NGC 6791. 

Recently, IBrogaard et al.l (1201 ll 120121) analyzed light 
curves of V18 and V20, two of the previously known eclips¬ 
ing binary systems in NGC 6791 with orbital periods_of 
18.8 days and 14.5 days, respectively dRucinski et al.|[l996l) . 


Their measurements of masses and radii of the individual bi¬ 
nary components are shown by open circles in the top left 
panel of Figure[9] Their reported ltr errors are 0.3% in mass 
and 0.6%-0.9% in radius, but 3cr measurement uncertainties 
in masses and radii are displayed as in lBrogaard et al.1 d2012i) . 
They also provided spectroscopic T e g- for the individual com¬ 
ponents, which are shown by filled circles in the_middle left 
panel with ±1 a error bars. In lBrogaard et al] (120121) . 7)n for 
the secondary of V20 was not directly determined from the 
spectra, but estimated assuming the same distance as for the 
primary. In the bottom left panel, filled circles show luminosi¬ 
ties computed from their measured radii and spectroscopic 
T e ff measurements. Errors in luminosity are those propagated 
from a 3cr error in radius and a la error in T e ff. 

The solid lines in the left panels of Figure [9] are theoreti¬ 
cal YREC models at 8.5 Gyr, 9.5 Gyr, and 10.5 Gyr, respec¬ 
tively, at the mean spectroscopic metallicity of NGC 6791 
([Fe/H]= ±0.37 ± 0.07). The two dashed lines additionally 
show 9.5 Gyr models with ±1 a changes in the adopted metal¬ 
licity. Models are shown for the MS only, since individual 
members of V18 and V20 are on the MS (Figure 0. As shown 
in the top left panel, the observed mass-radius relation can 
be matched satisfactorily with 9.5 or 10.5 Gyr models, when 
a 3d allowance is made in the comparison with the binary 
data. If one adopts a model that matches the mass-radius re¬ 
lation of V20 components (9.5 Gyr at [Fe/H]= ±0.37), both 
V18p and VI8s lie systematically above the model line, in¬ 
dicating that our model slightly underestimated stellar radii. 
The observed offset resembles earlier findings that standard 
one-dimensional stellar interior models generally predict ~ 
3% sm aller radii than actually measured for low-mass stars 
dSpada et al.l 120131 and references therein), although the bi¬ 
nary components in NGC 6791 are more massive than these 
stars (~ 0.2 M o -0.8 M 0 ). 

In the middle left panel of Figure |9j open circles indi- 
cate 7 d f inferred from the IRFM relations in C as a g ran d e et all 
(120101) . Since the IRFM relations in V — Ic are defined only 
up to [Fe/H]= ±0.3, we computed IRFM-based 7' c n at [Fe/H]= 
±0.37 from B - V colors (see iBrogaard et al.ll201 ll) for indi¬ 
vidual components in the binary systems. We took errors of 
0.02 mag in B — V, and corrected for the foreground extinc¬ 
tion using E(B — V) = 0.11 (see below). Except for the V20s, 
IRFM estimates are ~ 1 30-190 K cooler than the sp ectro- 
scopic values reported by Brogaard et al.l d201 ll I2012L filled 
circles), but lower temperatures show better agreement with 
our models. As shown in the bottom left pane l of Figure |9l 
larger radii and higher T e ff in IBrogaard et all (1201 ll 120121) 
yield larger luminosities of V18 components than our model 
predictions. 

In fact, the differences between the filled and open circles 
in Figure [9] essentially show a range of solutions available 
for eclipsing binary systems from independent 7)n estimates, 
since each of the methods is subject to systematic uncertain¬ 
ties in the parameter estimation. Photometric temperatures 
from IRFM depend on the foreground extinction correction, 
but r eff for both V18 components is still ~ 70 K cooler than 
spectroscopic estimates, even if we assumed the same red¬ 
dening [E(B - V) = 0.16] as in IBrogaard et al.l (120111) . The 
IRFM relations are also weakly constrained in this regime, 
where true empirical constraints would be welcome and are 
needed for cool, metal-rich stars. Spectroscopic methods pro¬ 
vide an independent check on photometric temperatures, but 
there could be systematic errors unaccounted for when sep¬ 
arating individual components of unresolved binary systems 
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Fig. 8.— CMDs of NGC 6791. The BVIc data are from IStetson et all 120031 ). and At-band data in the right panel are originally from ICarnev et alj (2005) 
after transforming their CIT K into K s magnitudes in the 2MASS system. Black plus points are stars remaining after the photometric filtering, which trace a 
single-star sequence of the cluster with little influence from cluster binaries and background stars. Open circles are cluster sequences derived from these points. 
Locations of individual members of VI8 and V20 eclipsing binary systems are marked with star symbols. The dotted diagonal line on the right panel represents 
a completeness limit in K s . 


from their spectra. Given the above limitations in each of 
these approaches, we conclude that our models are in reason¬ 
able agreement with the observed binary data in NGC 6791. 

The right panels in Figure [9] show the same binary data 
as in the left panels, but display a set of YREC models at 
8.5, 9.5, and 10.5 Gyr assuming AT/AZ = 2.0 (G. New- 
sham & D. M. Terndrup 2007, private communication; see 
also lNewsham & Terndruril2007l) . The corresponding helium 
mass fraction is Y = 0.32 at [Fe/H]= +0.37 (Z = 0.038), as 
opposed to Y = 0.29 in our base models with AT/AZ = 1.2. 
Comparisons with the binary data show that helium-enhanced 
isochrones better match data with spectroscopic temperatures 
(filled circles in the middle and bottom panels). However, 
AT/AZ = 2.0 is significantly steeper than our value con- 
strain ed fr om the primordial and the solar helium abundances 
(see § 12. lb . The amount of helium required is also larger than 
that of iBrogaard et all (120 121 T = 0.30 ± 0.01) based on an 
impro ved set of Victoria-Regina models (VandenBerg et al. 
120141) . Fortunately, the MS-fitting method, and thus our cal¬ 
ibration, is re lativ ely insensitive to the adopted AT/AZ as 
discussed in § 15.61 

Although there still remain unresolved issues about stellar 
radii, similar quality fits were obta i ned using other indepen- 
dent models (IBrogaard et all 1201 ll 1201 2t IVandenBerg et all 
120141) . The comparison with the binary data suggests that 
YREC interior models are in satisfactory agreement with ob¬ 
servations. In the next section, we apply color- T e ^ relations 
and the proposed empirical color-T e ff corrections to YREC 
models, and proceed to evaluate the accuracy of our color- 
calibrated YREC models by inspecting the internal consis¬ 


tency of MS-fitting results on CMDs with different color in¬ 
dices. The goal of this exercise is to verify that empirical 
corrections do not depend, or depend only weakly, on metal- 
licity, and that the corrections derived using the Hyades and 
Praesepe are valid for NGC 6791. 

5. TEST OF MODELS FROM MS FITTING 

In this section, we describe our effort to test the accuracy 
of the calibrated isochrones using CMDs of NGC 6791 in dif¬ 
ferent colors. We begin our discussion by inspecting the con¬ 
sistency of isochrone fitting in narrow color ranges, where a 
pair of cluster parameters (reddening versus metallicity and 
distance modulus versus metal licity ) can be constrained with 
a little dependence on others (§ 15. Ik We also use IRFM color- 
T e ff relations (§ 15.21 ), independent cluster photometry (§ 15.31 ), 
and cl uster photometry and models in the SDSS filter sys¬ 
tem (§ 15.41 ) to check the accuracy of our models and the input 
cluster photometry. We present a global search f or the cl uster 
parameters and robust error estimates in Sections [53] and [5T6J 
respectively. 

5.1. Preliminary Fitting with Calibrated Isochrones 

The top panel in Figure [TO] shows E(B -V) estimates as 
determined from a set of 9.5 Gyr old calibrated isochrones, 
spanning a range of input metallicity from +0.1 to +0.6 dex. 
Following the procedure described in Paper IV, we derived a 
color excess on each of BV, Vic , and VK S CMDs, by shift¬ 
ing a model along a reddening vector at a given distance. We 
restricted isochrone fits to 17.8 < V < 18.5, making our esti¬ 
mates of color excess insensitive to both age and distance. As 
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Fig. 9.— Comparisons of theoretical models with observed mas s-radiu s (top), mass-temperature (middle), and mass-luminosity (bottom) relations for VI8 and 
V20 eclipsing binary systems in NGC 6791 dBrogaard et all20TTll2012j) . Error bars are shown for 3cr errors in mass and radius, and lcr errors in T e ff. The error 
bars in luminosity include a lcr error in T e ff and a 3cr error in radius. Left: solid lines are 8.5, 9.5, and 10.5 Gyr old isochrones at the spectroscopic metallicity of 
NGC 6791 ([Fe/H]= +0.37), while dashed lines show 9 .5 Gyr o ld models with a d=l<+ change in metallicity (cr = 0.07 dex). In the middle and bo ttom panels, filled 
circles are measurements based on spectroscopic T e ff in IBrogaard et alj d2012D , while open circles show the case when T e ff is estimated using the lCasagrande et aH 
<201 OP IRFM relation at [M/H]= +0.37 and E(B — V) = 0.11. Right: same as in the left panels, but showing helium-enhanced models with AY/AZ = 2.0 (see 
text). 


described in § 12.11 we adopted extinction laws with an explicit 
dependence on colors, which produce slightly redder colors 
for MS stars in NGC 6791 than those without color terms: 
0.008, 0.005, 0.008, and 0.001 mag in B - V, V - I c , V - K s , 
and J - K s , respectively, when E(B - V) = 0.1. The E(B - V) 
estimates displayed in Figure[l0]are those for zero-color stars. 

We computed a distance using weighted median statistics, 
restricting model fits to 1.00 < B — V < 1.10 and 1.05 < 
V -I c < 1.15 (4900 K < Teff < 5200 K). This part of the 
MS has a relatively shallow slope on CMDs, where a distance 
can be constrained with little dependence on reddening. We 
used a mean distance modulus obtained from the BV and Vic 
CMDs when estimating a color excess on the VK S CMD. In 
this way, we iteratively solved for both color excess and dis¬ 
tance of the cluster. The results are shown as triangles, boxed 
points, and cross marks for BV, Vic, and VK S CMDs, respec¬ 
tively. The solid, dotted, and dashed lines are a linear fit to 
MS-fitting results on each CMD, as a good approximation of 
the observed trend in the above metallicity range. 

In the top panel of Figure [TO] the line from the BV CMD 
has a slope that is steeper than those from the other two 
color indices. In this preliminary fitting exercise, we used 
differential sensitivities of colors on metallicity (B - V ver¬ 
sus V — Ic and B - V versus V - K s ) to constrain the clus¬ 
ter’s metallicity and foreground reddening (see Papers III 
and IV). The line from the BV CMD crosses the Vic CMD 
line at [M/H]= +0.459 ± 0.012, and the reddening at this 
metallicity corresponds to E{B — V) = 0.095. On the other 


hand, the intersection between the BV and VK S lines has 
[M/H]= +0.369 ± 0.009 and E(B - V) = 0.128. A weighted 
mean photometric metallicity from these two color combi¬ 
nations yields [M/H]= +0.40 ± 0.05, which is in excellent 
agreement with the spectroscopic metallicity of the cluster 
([Fe/H]= +0.37 ± 0.07). The error represents half of the dif¬ 
ference between the two photometric estimates, while an error 
propagated from fitting a model to a CMD is about six times 
smaller than this value. At the mean photometric metallic¬ 
ity, the reddening of the cluster is E(B — V) = 0.115 ± 0.010 
(scatter from the three CMDs) ±0.009 (propagated from the 
photometric metallicity error). 

Our photometric solution for reddening can also be con¬ 
strained by adopting the accurate cluster metallicity in the 
literature. At the mean spectroscopic abundance ([Fe/H ]= 
+0.37 ± 0.07), as shown by vertical lines in Figure IIOFI the 
mean reddening from BV, Vic, and VK S CMDs becomes 
E(B - V) = 0.121 ± 0.011 (scatter from the three CMDs) 
±0.014 (propagated from the metallicity error). The disper¬ 
sion in E(B -V) from the three CMDs may reflect small sys¬ 
tematic errors unaccounted for in the current isochrone cali¬ 
bration. However, other sources of systematic errors are also 
likely. For example, ICarnev et all ( 2005 ) made a comparison 

8 Since a number of spectroscopic studies adopted color-7 e ff relations to 
calculate temperatures of stars, their abundances are correlate d w ith a fore¬ 
ground reddening of the cluster adopted in each study. In Figure lTol we simply 
use vertical lines to represent a range of spectroscopic metallicity determina¬ 
tions, neglecting a correlation between [Fe/H] and E(B — V). 
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Calibrated Models 



[M/H] 


Models with Lejeune et al. Colors 



Fig. 10.— Top: color-excess estimates for NGC 6791 from CMDs with 
B - V (triangles), V — Ic (boxes), and V — K s (crosses) using 9.5 Gyr old 
calibrated isochrones at various metallicity bins. Lines are a linear fit to each 
set of data points. Bottom: same as in the top p anel, but using models with 
original color-7^ transformations in lLeieune et aljfi 9971119981) . The vertical 
lines represent the spectroscopic metallicity of the cluster and its ±ler errors. 


of their photometry with bright 2MASS stars, and found that 
their transformed At,-band photometry is 0.014 ±0.005 mag 
brighter than 2MASS. If we allow this change to V - K s col¬ 
ors, E(B — V) from the VK S CMD would become smaller by 
A E(B — V) = 0.006, reducing the scatter in the estimates of 
color excess from the three CMDs at the spectroscopic metal¬ 
licity of the cluster. 

Our E(B — V) is found within a broad range of previous es¬ 
timates in the literature (see §Q}. Reassuringly, it is also com¬ 
parable to or slightly lower than the integrated value [E(B — 
V) = 0.1331 from the dust emission map (Schlegel 'et al.ll 19981: 
ISchlaflv et all l2010h . The cluster is located high above the 
Galactic plane (|z| ~ 750 pc), which is at least five times the 
scale height of the dust layer in the Galactic disk. There¬ 
fore, most of the dust should be found between the Sun and 
the cluster, and the integrated line-of-sight reddening should 
serve as a good approximation for the cluster reddening. 

We also attempted to estimate E(B -V) from the JK S CMD, 
but our values were always too high compared to those ob¬ 
tained using other color indices: at [M/H]~ +0.4, the JK S 
CMD yields E(B - V) ~ 0.24 ± 0.01. Since E(J - K S )/E(B - 
V ) = 0.56, A E(B — V) = 0.12 mag requires a ~ 0.07 mag 
change in the zero point of J - K s colors. While a transforma¬ 
tion from the UKIRT JIIK system to the CIT or 2MASS sys¬ 
tem is not trivial, Carney et al . ( 200 5b found that their J - K s 
photometry is only 0.007 ±0.005 mag redder than 2MASS. 
The comparison with IRFM (Figure |6]> shows that our model 
J - K s colors are too blue by 0.02-0.03 mag, but the differ¬ 
ences are still too small to explain the large color excess on 


1.0 ^ B-V ^ 1.1 



Fig. 11.— Distance estimates for NGC 6791 from CMDs with B-V (tri- 
angles) and V — Ic (boxes) colors in various metallicity bins. The 9.5 Gyr 
old isochrones with empirical color corrections were used in two different 
color bins: 1.0 <B — V < 1.1 (top) and 1.2 < B — V < 1.4 (bottom), and 
corresponding color ranges in V - Ic- The foreground reddening was fixed at 
£(£-V) = 0.11. 

the JK S CMD. This may indicate a problem either in the IR 
color-r e ff relations at high metallicity or in the zero point of 
the photometry. We leave the discussion on the calibration of 
IR colors to the next paper of this series. 

In contrast, models without empirical color corrections re¬ 
sult in a significant inconsistency in the derived reddening as 
shown in the bottom panel of Figure [lO] The E(B - V) es¬ 
timates from the BV CMD are systematically higher, while 
those from the Vic and V K s CMDs are lower than those ob¬ 
tained using the calibrated isochrones (top panel). Our empir¬ 
ical corrections do not include metallicity terms, and therefore 
the slopes of the lines in the bottom panel are similar to those 
in the top panel. On the other hand, the intercepts of these 
lines are modified by our color-r e ff corrections, leading to an 
improved internal consistency of E(B - V) estimates at the 
spectroscopic metallicity of the cluster (top panel). 

Similarly, the top panel in FigurefTTIshows a distance mod¬ 
ulus as a function of metal abundance in the model. We fixed 
the foreground reddening at E(B - V) = 0.11 to be consistent 
with the above result and assumed an age of 9.5 Gyr. The dis¬ 
tance modulus was derived using stars with 1.00 < B - V < 
1.10 and 1.05 < V - Ic < 1.15 in the BV and Vic CMDs, 
respectively. Triangles show distance moduli from the BV 
CMD, while boxed points are those from the Vic CMD. We 
excluded V — K s because its CMD only covers the vertical 
part of the MS. Solid lines are a linear fit to each set of the 
data points. At higher metallicities, models become brighter, 
and a greater distance is derived from MS fitting. However, 
the BV CMD is more sensitive to a metallicity change than the 
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Vic CMD, and consequently has a steeper slope in Figure fill 

At the photometric metallicity of the cluster ([M/H]= +0.40; 
see Figure [Toll, isochrone fitting in the top panel of Figure fill 
yields (m — M)o= 13.074±0.005 and 13.006±0.005 from the 
BV and VI C CMDs, respectively, where the errors represent a 
scatter of data points from the best-fitting isochrone. Within 
the fitting errors alone, distances from these two color indices 
are inconsistent with each other. As noted above, this differ¬ 
ence can originate from small systematic errors in our models, 
but it can also be easily ca used by various systematic errors in¬ 
volved in the MS fitting (§ 15.61 ). The two lines cross at a lower 
metallicity than the spectroscopic metallicity of the cluster, 
indic ating that the photometric metallicity is [M/H]« +0.26. 
In § 15.51 we seek a global solution of the cluster’s param¬ 
eters from all available features on CMDs, where the dif¬ 
ference from a model in these color indices is minimized. 
At the spectroscopic metallicity of the cluster, the average 
distance modulus from the B — V and V — Ic color indices 
is (in — M) o = 13.002±0.042 (from the difference between 
the two CMDs) ±0.067 (propagated from the metallicity er¬ 
ror). Our distance estimate is within the range found in pre¬ 
vious studies in the literature, 12.9 < (m — M)n < 13.1 (e.g., 
Salaris_etalj|2004yCarne y et al .1120051 : 1 An thonv-Twarog et all 

20071: iBrogaard et alJl201 1L~ among others). 

However, the mutual agreement of our models in B — V and 
V - Ic breaks down when we go to cooler stars. The bot¬ 
tom panel in Figure |TT| shows the case when fits were made 
using redder stars in NGC 6791 (1.20 < B -V < 1.40 and 
1.25 < V - I c < 1.65, or 4300 K < 7; ff < 4750 K) than those 
shown in the top panel. A strong inconsistency between B — V 
and V - Ic colors is seen even after the process of empiri¬ 
cal color correction. Distance moduli from the BV CMD are 
systematically larger by ~ 0.1 mag than those obtained from 
bluer stars, while distance estimates in V — Ic remain almost 
unaffected. This implies that systematic errors in B — V col¬ 
ors, either in the photometry or in the model in the above tem¬ 
perature range, are significantly larger than in V — Ic- This is 
unfortunate, since fitting along the MS, in the presence of such 
systematic offsets, could reduce the accuracy in our distance 
and metallicity estimates considerably. 

The left panels in Figure fl2l show color residuals from the 
best-fitting isochrones for individual stars in B - V (top pan¬ 
els) and V - Ic (bottom panels). We used 9.5 Gyr old models 
at the spectroscopic metallicity of the cluster ([Fe/H]= +0.37). 
Colors were corrected for the reddening with E(B-V) = 0.11, 
where we adopted a distance modulus as derived using stars in 
1.0 < fi - V < 1.1 on each CMD: (m - M) 0 = 13.037 in BV and 
(m - M )o = 12.983 in the Vic CMD. Only those data points 
that were left from the photometric filtering routine (see Fig¬ 
ure [Sj are shown. 

As seen in the top left panel of Figure [12] there is a signif¬ 
icant residual trend in B - V at 19.5 < V < 21.5. Large cir¬ 
cles represent moving-averaged differences in bins of AT = 
0.5 mag, which are tabulated in Table [4] The largest color 
deviation (~ 0.03 mag) is seen at V ~ 20.75 mag or at 
B - V ~ 1.3, in the sense that our best-fitting model is red¬ 
der than the observed MS. Given that the slope of the MS is 
about 4.5 in the BV CMD, an average offset of ~ 0.02 mag 
in B — V is consistent with the ~ 0.1 mag larger distance 
moduli obtained using redder stars in Figure [TT] The ob¬ 
served offsets persisted even if we adopted different metal- 
licities (+0.3 < [M/H] < +0.5) of models. On the other 
hand, the V — Ic color differences in the bottom left panel are 
found within A(T — Ic) < 0.015 mag, down to V ~ 22 mag 


(V - Ic ~ 2.0). Even at the maximum displacement in B - V, 
the match to the data in V - Ic is excellent. That the model 
matches the observed MS in the Vic CMD relatively well sug¬ 
gests that the systematic residuals in B — V were probably not 
caused by errors in the stellar interior models. 

In order to check that the large residual trend in B - V was 
not caused by the adopted reddening of the cluster, we show 
color residuals for an alternative reddening [E(B - V) = 0.14] 
in the right panels of Figure [12] As in the left panels, we 
derived and adopted distances from individual CMDs using 
stars in 1.0 < B — V < 1.1: (m - M) o = 13.096 in BV and 
( m - M) o = 13.092 in the Vic CMD. As seen in these panels, 
the large B - V color residuals do not disappear, but color 
residuals in V — Ic become significantly large. The V — Ic 
color index reacts more sensitively to E(B — V), suggesting 
that E(B - V) = 0.14 is too high. 

There are a few remaining possibilities that could have 
made the observed MS in B — V deviate from the model. 
Firstly, the systematic color offset could originate from in¬ 
accurate color- T e ff relations in the model. The higher metal¬ 
licity causes stars to have greater line blanketing at shorter 
wavelengths. The effect of line blanketing would be en- 
hanced_at lower t emper atures, so we can hypothesize that 
the [Lejeune et al . ( 1199 7!) relations failed to reproduce fluxes 
for cool MS stars in NGC 6791 because of the difficulty in 
modeling complicated spectral energy distributions in shorter 
wavelength passbands, such as in B. The effect of line blan¬ 
keting would be less in V and certainly less in Ic, so V - Ic 
colors would be less affected. Second, it is possible that the 
problem lies with systematic errors in the cluster photome¬ 
try. IStetson et al] (120031) claimed that zero-point errors of 
their photometry are at least 0.0014 mag in V and Ic, and 
0.0025 mag in /(; the corresponding zero-point error in B — V 
is at a minimum of 0.003 mag. In addition, there could ex¬ 
ist color-scale errors in the photometry. As mentioned above, 
the expected amount of the color-scale error is ~ 0.02 mag 
in B - V. Below we proceed to look for effects that could 
produce the discrepancy by using alternative color- T e ff trans¬ 
formations and independent cluster photometry. 

5.2. IRFM Color-T e ff Relations 

We utilized color-T e ff relations in ICasagrande et all (120101) 
to test whether the B - V color offset is found independently 
of models employed. Figure[l3]is the same as FigureQT] but 
shows IRFM-based slopes and intercepts. We used the same 
base isochrones, but after transforming model T e ff into B — V 
and V — Ic using the lCasagrande et akl (12010 ) relations. Given 
that their color-T e ff relations are mostly empirical with little 
dependence on models, an extra step of empirical calibration 
as described in this paper is not needed. The same set of stars 
as in Figure[TT]was used in the top (1.00 <B -V < 1.10) and 
the bottom (1.20 < B - V < 1.40) panels. 

In the top panel, the IRFM method differs from what our 
models predict in that the metallicity dependence and the 
zero points are slightly different. The lines have slopes of 
1.16 ± 0.02 mag dex -1 for B - V and 0.62 ± 0.04 mag dex -1 
for V — Ic, while those in our models (top panel in Figure ITTli 
have 1.20 ±0.01 mag dex -1 and 0.72 ±0.01 mag dex -1 , re¬ 
spectively. As a consequence, the IRFM-based photometric 
solution to metallicity is somewhat different from what we 
obtained above, and yields [Fe/H]=+0.30. Since it is within 
a lcr bound of the spectroscopic metallicity of the cluster, we 
consider that there are no alarming signs of mismatch with the 
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Fig. 12.— The B — V and V — 1q color residuals of stars in NGC 6791 in the sense of observed colors minus calibrated models at [M/H]= +0.37. The left and 
right panels show comparisons with E(B — V) = 0.11 and 0.14, respectively, at a distance derived on each CMD. Only those points left from the photometric 
filtering are shown. Large circles are average color differences in bins of AV = 0.5 mag. Their error bars are smaller than the size of a symbol. 




TABLE4 

Average Color Differences from Calibrated Isochrones in NGC 6791 


V 

E(B - V) 

= 0.11 

E(B - V ) 

= 0.14 

(mag) 

A JB - V) 

All/ - l c ) 

A (B-V) 

A(V - lc) 

17.75 

0.0047 ± 0.0007 

-0.0013 + 0.0009 

-0.0313 + 0.0011 

-0.0467 ± 0.0011 

18.25 

0.0080 + 0.0005 

-0.0033 + 0.0007 

-0.0069 + 0.0006 

-0.0256 + 0.0007 

18.75 

0.0059 + 0.0005 

-0.0060 + 0.0008 

0.0002 ± 0.0005 

-0.0151+0.0008 

19.25 

-0.0007 + 0.0006 

-0.0025 + 0.0008 

-0.0008 + 0.0006 

-0.0033 + 0.0008 

19.75 

-0.0154 + 0.0007 

0.0124 + 0.0009 

-0.0107 + 0.0007 

0.0216 + 0.0010 

20.25 

-0.0256 + 0.0009 

0.0012 + 0.0010 

-0.0192 + 0.0009 

0.0346 + 0.0009 

20.75 

-0.0283 + 0.0012 

-0.0053 + 0.0011 

-0.0235 + 0.0012 

0.0278 + 0.0011 

21.25 

-0.0226 + 0.0019 

-0.0027 + 0.0014 

-0.0239 + 0.0019 

0.0373 + 0.0014 

21.75 

-0.0035 + 0.0031 

-0.0046 + 0.0019 

-0.0147 + 0.0031 

0.0320 + 0.0016 

22.25 

0.0032 + 0.0048 

0.0251 + 0.0027 

-0.0169 + 0.0048 

0.0496 + 0.0029 


NOTE. — The sense is observed minus model colors. Errors represent a standard error in the mean color 
difference. 


data. 

As shown in the bottom panel of Figure [T3j however, a 
significant discrepancy in distance is seen between BV and 
Vic CMDs when fits were limited to cooler stars (1.20 < 
B - V < 1.40 or 4300 K < 7 dr < 4750 K). The observed be¬ 
havior is similar to those seen in FigureQT]with our calibrated 
isochrones. In other words, the two independent color-r e ff 
relations show that a distance modulus from the BV CMD 
strongly depends on the color range chosen for MS fitting, 
but the distance from the Vic CMD is certainly less affected 
by the choice. 

Figure [14] shows color differences from a 9.5 Gyr old 
model at [M/H]= +0.37 with E(B - V) = 0.11, after trans¬ 
fo rming model 7j,ff into c olors using the IRFM relations 
in iCasagrande et alJ d2010l) . The comparison is limited to 
V ~ 20.5 mag due to the T e ff limit set in the IRFM work. 
We used a distance modulus [{m — M )o = 13.010] that best 
matches the model with the IRFM relation on the BV CMD 


in 1.00 <B-V < 1.10. In the right panel, we also show the 
case when an alternative value of reddening [E(B - V) = 0.14] 
was adopted at its best-fitting distance [(m - M)o = 13.048]. 
As seen in Figure [14] the systematic residual trend from the 
IRFM-based models also shows that the B - V photometry is 
too blue for fainter stars (V > 19.5 mag), independently of the 
cluster reddening adopted. Interestingly, the systematic trend 
is even stronger than those seen in Figure[l2]based on our cal¬ 
ibrated models. This is not surprising, as the empirical IRFM 
relations are not as well defined for extreme metallicity and 
temperature. 

A strong internal inconsistency of IRFM colors for cool 
dwarfs can also be seen without relying on YREC interior 
models. Figure[T5lcompares 7' c n estimates from B — V,V — K s , 
and J - K s with those from V — Ic based solely on IRFM re¬ 
lations. Only those remaining from the photometric filter¬ 
ing (black cross points in Figure [8J are shown. Observed 
colors were corrected for the cluster’s foreground redden- 
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IRFM: 1.0 ^ B-V ^ 1.1 



[ M/H ] 

IRFM: 1.2 ^ B-V ^ 1.4 



[M/H] 


Fig. 13.— Sa me as in Fig ure 1 11 1 but distance moduli estimated using our 
base models with ICasagrande et alj (201 CD IRFM relations. 

ing, b efore T^f were estimated using the Casa grande et alJ 
(2010) relations. We assumed E(B - V) = 0.11 in the left 
panels and E(B -V) = 0 .14 in the right panels. Because the 
ICasagrande et a 1.1 (1201 Ol) relations were defined up to [Fe/H]= 
+0.3 in V - /(-. we assumed [Fe/H]= +0.3 in the above com¬ 
parisons. However, metallicity effects are relatively small in 
the differential T e g comparisons. For V — K s and J - K s , we 
only included stars with V < 18.5 mag to avoid a bias below 
the detection limit in K (dashed line in Figure [8]). 

Figure [15] clearly shows that there exists a strong system¬ 
atic trend in the temperature difference (Ar e ff) between B-V 
and V — Ic colors. Temperatures inferred from B-V are on 
average 80 K higher than those from V — l c at T e ff < 4800 K, 
while temperatures from B-V and V - Ic are found on nearly 
the same scale (< 10 K) at T e ff > 5000 K. The internal incon¬ 
sistency of 7cn persists even at the spectroscopic metallicity 
of the cluster. From the metallicity sensitivity of the IRFM 
relations, we estimated that a ~ 0.1 dex increase in [Fe/H] 
results in a ~ 30 K increase in temperature from both B — V 
and V — Iq colors, but the change in metallicity does not make 
any appreciable change in the systematic 7c it difference. The 
higher temperatures in B-V than in V — Ic at 7' c ff < 4800 K 
imply that B — V colors are too blue at a given V - Ic color, 
and are consistent with systematically larger distance moduli 
from B — V in Figure[TT] Unless both our and IRFM color-T e ff 
relations are in error, the internally inconsistent T s s suggests 
that the problem lies in the cluster photometry of cool MS 
stars. 

5.3. Other Cluster Photometry 

We compared Stetson’s c lu ster photometry with those in 
iKaluznv & Rucinskil (119951) . iKaluznv & Rucinskil obtained 
deep UBV photometry of stars in the central part of the cluster 


E(B-V) = 0.1 1 




Fig. 14.— Same as in Figure fl2l but showing co lor res iduals in B-V from 
a [M/H]= +0.37 model with the lCasagrande et all (2010 5 IRFM relation. The 
left and the right panels show comparisons at E(B — V) = 0.11 and E(B — 
V ) = 0.14, respectively. A comparison in V — Ic is not shown because of the 
metallicity limit in the IRFM relation. 

using the KPNO 2.1 m telescope, and took shallower UBVIc 
data over the entire cluster field using the KPNO 0.9 m tele¬ 
scope. Below, the former data set is referred to as Table 1, 
and the latter as Table 2, following the file numbering con¬ 
ventions in the VizieR archivesQfrom which we downloaded 
their photometric tables. We used a matched coordinate list in 
their Table 2 and relatively bright stars in lStetson et all (120031) 
to derive our own astrometric plate solutions and to transform 
pixel positions into celestial coordinates for their data sets. 

Figure fl6l shows magnitude and color differe nces between 
Table 1 in IKaluznv & Rucinskil (119951) and iStetson et all 
(120031) . in the sense of the former minus the latter. 
Similarly, Figure [IT] shows differences from Table 2 in 
IKaluznv & Rucinski The blue diamond point is a median 
difference in 0.5 mag intervals in V, 0.05 mag in B-V, and 
0.06 mag in V - Ic- Error bars represent a median of all ab¬ 
solute deviations (MAD) multiplied by 1.483 to make the es¬ 
timator consistent with the standard deviation for a normal 
distribution. The scatter of data points is large, but a stan- 

s http://cdsarc.u-strasbg.fr/ 





















































Distances from Main-Sequence Fitting. V 


17 


E(B-V) = 0.1 1, [Fe/H] = + 0.3 


200 

0 




-200 - 


400 

200 

0 

-200 

-400 

800 

400 

0 

-400 F- 
-800 r 


T„ f (B-V) - T.„ (V-I c ) 


T e „ (V-K.) - T eff (V-I c ) 


T e ff (4 K s ) T ef( (V . I c ") 



4500 


5000 

■ (V - l c ) 


200 

0 


E(B—V)=0.1 4, [Fe/H] = +0.3 


*"X ! '% M\ 


-200 - 


400 

200 

0 

-200 

-400 

800 

400 

0 


-400 r 


-800 r 


5500 


Teff (B-V) - T eff (v-l c ) 



4500 


5000 

■ (V - Ic) 


5500 


FlO. 15.— Differences in T e ff as inferred from the lCasagrande etaD <2010 1IRFM relations in B — V, V — Ic, V — K s , and J — K s . Only those selected as single 
cluster members in NGC 6791 are shown. Observed colors were corrected for the cluster foreground reddening with E(B — V) = 0.11 and 0.14 in the left and 
right panels, respectively. Because the IRFM relation in V — Ic is defined at [Fe/H]< +0.3, temperatures were estimated at [M/H]= +0.3 in all color indices. 


dard error in the mean difference (a/y/N) is typically smaller 
than the size of a diamond symbol. In these c omp arisons, 
we imposed the same selection criteria as in IStetson et al.L 
who have also done star-by-star comparisons with photom¬ 
etry in lKaluznv & Rucinskil These were based on sharp, 
separation indice s, photometric errors, as well as a photo¬ 
metric quality flag in iKaluznv & Rucinski, Instead of select- 
ing stars with good quality flags in both IKaluznv & Rucinskil 
tables, however, we only used quality flags from each ta¬ 
ble when comparing it with IStetson et al.i photometry. This 
seemed to make more stars available in the photometric com¬ 
parisons, especially fainter stars at V > 20, which were not 
included in the original comparisons in IStetson et all 
Figures [16] and |T7] c learly show t hat there are systematic 
errors present in either IStetson et all or IKaluznv & Rucinskil 
photometry. For stars near the cluster’s turn-off (V ~ 17.5 or 
B-V ~ 0,9), differences in V are of the order of 0.02 mag. 
Besides photometric zero-point errors. Figures [16] and 1771 
demonstrate that differences in photometry depend on the 
magnitude or color of a star. Near the cluster’s turn off, 
B-V colors agree with each other among the three data sets. 
However, the difference increases toward fainter magnitudes 
and redder colors for both IKaluznv & Rucinskil tables, which 
amounts to 0.02 mag at B - V ~ 1.3. Interestingly, this offset 
can explain the difference we found in Figure|TT]between our 
calibrated isochrones and the IStetson et all photometry. Al¬ 
though it is difficult to make a fair judgment about which pho¬ 
tometry is correct based on the photometric comparison alone. 


one might imagine that a color-scale error in the Stetson’s 
photometry has not been completely removed in the newer 
data reduction (IStetsonll2005l) . If Stetson’s B-V colors are 
systematically bluer at fainter magnitudes, higher IRFM tem¬ 
peratures in B — V than those from V — Ic (Figure fl4l) can be 
naturally explai ned. The V - Ic co lors also show large d iffer- 
ences between IStetson et all and IKaluznv & Rucinski (bot¬ 
tom panel in Figure fTTI). However, our models show no signs 
of large systematic color residuals in V - Ic (left panel in Fig¬ 
ure Eh and distances derived using bluer stars are consistent 
with those from redder ones (Figure [TTh 

5.4. Co/or-r e ff Relations in SDSS Colors 

We further checked the accuracy of model colors in the 
SDSS system. Figure[l8]is the same as FigureQTl but show¬ 
ing metallicity sensitivities in SDSS colors. We took griz pho¬ 
tometry of the cluster from I An et al . (200 81). As described in 
lAn et al.L the cluster photometry was extracted from the origi¬ 
nal SDSS imaging frames u sing th e DAOPHOT/ALLFRAME 
suite of programs (jStetsonlll987l 119941) . and was put on the 
natural ugriz system in SDSS. After applying photometric fil¬ 
tering on gr, gi and gz CMDs, we employed isochrones in 
lAn et all ( 20091). see Table H} to derive a distance modulus on 
each CMD using the r-band as a luminosity index. The inte¬ 
rior models are the same as those used in the current study, 
but the colors are those corrected based on observations of 
M67, on top of theoret ical colors computed using MARCS 
(iGustafsson et al.ll2008l) . We assumed the same cluster age 
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KR95 Table 1 - Stetson 
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Fig. 16.— Co mparisons of the revis ed Stetson et al. photometry with val- 
ues in Table 1 inlKa luznv & Rucinski (1995). The difference is in the sense 
IKaluznv & Rucinskil values minus those in lStetson et all The blue diamond 
points are median differences, and the error bars indicate a median absolute 
deviation of the difference multiplied by 1.483. A standard error in the mean 
difference is typically smaller than the size of a diamond symbol. The red 
dashed line represents a zero difference. 


KR95 Table 2 - Stetson 



V—l c (mag) 

Fig. 17.— Same as in Figure lT6l but showing comparisons with photometry 
in Table 2 in lKaluznv & Rucinskil <19951) . 
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[M/H] 
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Fig. 18.— Same as in FigureQT] but showing MS-fitting results from the 
application of calibrated isochrones in the SDSS griz colors JAn et al 1200914) - 


and reddening as in Figure [IT] 

Figure fl8l shows results when our fits were limited to the 
same color ranges as in FigureQT] In the top panel, MS fitting 
included stars in 0.80 <g-r< 0.95, which is similar to 1.0 < 
B-V <1.1, and the corresponding color ranges in g — i and 
g — z. In the bottom panel, results are shown for stars in 1.05 < 
g — r <1 .30 or 1.2 < — V < 1.4. In both panels, we assumed 
E(B - V) = 0.11. Unlike in the Johnson-Cousins and 2MASS 
systems, differential sensitivities in griz colors are limited; in 
other words, the lines from different colors in Figure [18] have 
similar slopes, leading to a large uncertainty in the estimate of 
photometric metallicity (see lAn et al.ll2009~ai.l2013 ). 

Nevertheless, Figure QT] clearly shows that the distance es¬ 
timates from the three SDSS color indices are in good agree¬ 
ment with each other at the spectroscopic metallicity of the 
cluster. In addition, almost the same distance modulus of 
the cluster is obtained independently of color ranges chosen 
for the fits. At the spectroscopic metallicity of the cluster, 
the mean distance moduli from the three color indices are 
(m-M)o= 13.009±0.034and 12.994±0.050from the blue 
(top panel) and the red (bottom panel) color fitting ranges, 
respectively. The errors are dominated by the error in the 
spectroscopic metallicity, while individual fitting errors and 
the dispersion in distance among the three SDSS color in¬ 
dices are small (a < 0.01 mag). Within the errors, both dis¬ 
tance estimates are also in excellent agreement with our BVIc- 
based photometric solution from the blue part of the CMD 
[(m — M ) o = 13.00±_see top panel in Figure [TIT In FigureQT] 
the lStetson et ali ( 2003 ) photometry shows a maximum devia¬ 
tion from the theoretical model at B - V ss 1. 3 (T e p - ~ 4600 K), 
which corresponds to g — r ss 1.1 (iAn et aljl2013l ). However, 
such an offset is not seen on the gr CMD (see Figure 18 in 
[An et al.ll2009bh . 

To summarize, the good agreement between photometric 
and spectroscopic metallicity (top panels in Figures [TO] and 
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0 suggests that our isochrones in B - V and V — Ic are mu¬ 
tually consistent after our process of empirical color correc¬ 
tion. The distance and reddening we found are also close to 
values obtained in the literature. However, the mutual agree¬ 
ment breaks down for cooler stars (lower panel of Figure fill. 
Isochrones using empirical corrections are too red in B — V 
compared to the photometry of stars in NGC 6791. We also 
found a similar color offset using IRFM color-r e ff transforma¬ 
tions. The observed offset could be due to a color-scale error 
in the cluster photo metry, as we found that the B — V photome¬ 
try inl Stetson et al. ( 2003 ;! is systematically bluer for cool MS 
stars in comparison with alternative cluster photometry. Al¬ 
ternatively, the observed difference could be caused by differ¬ 
ences between Stetson’s and our compiled cluster photometry 
of the Hyades and Praesepe. It can also originate from errors 
in the models, but we obtained internally consistent fitting 
results in the SDSS colors using models with similar input 
physics and color calibration. To avoid the large color offset, 
we excluded B — V photometry of cool stars (B-V > 1.10) in 
the following estimation of the cluster parameters. 

5.5. Joint Determination of Cluster Parameters 

As discussed above, differential sensitivities of broadband 
colors can be used to constrain the reddening, distance, and 
metallicity of NGC 6791. In this section, we simultaneously 
solved for all these parameters and the age of the cluster, 
by directly comparing calibrated models with the observed 
cluster sequence (open circles in Figure [8j. The fit covered 
all parts of the CMDs including specific features that were 
used to const rain foreground reddening, metallicity, and dis¬ 
tance (§ ED , but excluded B-V photometry of cool stars 
[B — V > 1.10) as discussed above. In the next section, we 
provide a robust error analysis on our parameter estimates. 

To simultaneously determine the cluster’s metallicity, red¬ 
dening, distance, and age, we searched for the minimum y 2 
value as defined below: 

3 Nj 3 Nj ,,, _ \2 /y _ Y \2 

2 _ ^ 2 _ ^ \Vij ’'model, j) , vM j ^modeljV 

Xtot “ 2_^ 2_^ X-ij ~ 72 

j= 1 ,■= 1 j= 1 ,=1 L vii a x,ij J 

(i) 

where the subscript i indicates individual points of the ob¬ 
served sequence (open circles in Figure [8j on the jth CMD 
(each with B — V,V — Ic, and V - K s , respectively). Nj is the 
number of points included in the fit. V t j is a star’s V magni¬ 
tude, and Xjj is its color, while y mo deij and Ainodci./ represent 
the magnitude and color of a model, respectively. The oy,ij 
and (Tx.ij are photometric errors in V magnitude and color, 
respectively. 

We performed a grid search in [M/H], E(B - V), and age 
to find a minimum \ 2 ot . The average distance modulus was 
determined using the BV and Vic CMDs. In the distance 
estimation, we restricted model fits to B — V < 1.10, but in¬ 
cluded redder stars in the Vic CMD with V - Ic < 2.0. The 
V - Ic = 2.0 limit was set to avoid regions where our empir¬ 
ical corrections are poorly defined. We further restricted fits 
toV >18 mag in the distance estimation to exclude subgiant 
stars in the cluster. Adding the VK S CMD in the fitting pro¬ 
cess helped to better constrain the cluster’s age, metallicity, 
and reddening with its longer wavelength baseline (see the top 
panel in Figure flQt. although it was not used in the distance 
estimation. 

The top and bottom panels in Figure [19] show likelihood 
contours in [Fe/H] versus E(B — V) and [Fe/H] versus age, 



Fig. 19.— The Ax 2 distribution from MS fitting for NGC 6791. Likeli¬ 
hood contours are shown at Ax 2 0t = 4.7, 9.7, 16.3, and 23.5 (68.3%, 95.4%, 
99.73%, and 99.99% confidence levels for four degrees of freedom) relative 
to a minimum x 2 0 t- Error bars indicate ±1 a total errors in Table [6] at the 
minimum x 2 0 r The solid and dotted vertical lines represent the spectroscopic 
metallicity of the cluster and its ± la errors. 

respectively. The \ 2 surfaces were smoothed using a Gaus¬ 
sian kernel. Contours are shown at A\f ot = 4.7, 9.7, 16.3, and 
23.5 (68.3%, 95.4%, 99.73%, and 99.99% confidence levels 
for four degrees of freedom) relative to the minimum value of 
X 2 0t . The minimum Xtot va lue was large (xtot ~ 1000), sug¬ 
gesting that photometric errors were underestimated and/or 
that there were remaining systematic errors in the model cal¬ 
ibration. For a conservative error estimate, we normalized 
photometric errors by multiplying by a factor of 3.5 in all 
CMDs (see Paper III) to make x 2 lt « 89, which is approx¬ 
imately equal to the total number of degrees in the fitting 
{v = 96-4 = 92). 

At the global yf minimum, our solution yields 
[M/H]= +0.38, E(B - V) = 0.108, 9.5 Gyr, and a corre¬ 
sponding cluster distance modulus is (m - M) 0 = 13.00. 
Our best-fitting cluster parameters are in go od a greement 
with those from our earlier discussions in § ED Impor¬ 
tantly, Figure [T9] shows that our photometric metallicity 
estimate is in excellent agreement with the spectroscopic 
metallicity of the cluster. Our best-fitting age based on non¬ 
diffusion models (9.5 Gyr) could be ~ 1 Gyr too old (e.g., 
Proffitt & Mtndenberd [1991t ICastellani & degl’Innocentil 
199Pi: Ichahovcr et all 1200it) compared to the case with 
microscopic diffusion, although an accurate age scale is not 
the main concern of this paper. 

We adopted a conservative 1 er fitting error from the size of 
contours at Axf ot = 9.7, yielding cr([Fe/H]) = 0.015, a[E(B - 
y)] = 0.004, and 0.2 Gyr in age. For a la fitting error 
in distance modulus, we propagated these errors assuming 
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TABLE 5 

Summary of Cluster Parameters for NGC 6791 



[M/H] 

Age 

E(B-V) 

(m-M)o 

Solution 

(dex) 

(Gyr) 

(mag) 

(mag) 

Global x 2 Minimum 

0.38 + 0.06 

9.5 ±0.3 

0.108 + 0.014 

13.00 + 0.07 

Bias-corrected for Binaries 

0.42 ± 0.07 

9.5 ±0.3 

0.105 + 0.014 

13.04 + 0.08 


A(m — M)o/A[Fe/H] « 1 and A{m - M)o/AE(B — V) ~ 2 
(Paper III), along with a fitting error on each CMD. However, 
final errors are dominated by the sum of various systematic 
errors as described in the next section. Best-fitting cluster 
parameters from the global minimization and their total 
errors are listed in the first row of Table [5] 

Figure [20] displays CMDs with B — V, V — Ic, and V - K , 
indices, where the best-fitting isochrones are shown as solid 
red lines. Models with ±0.5 Gyr age differences are also 
displayed to show the quality of fits with varying ages. The 
dotted-dashed b lue fines are the same theo retical model, but 
with the original lLeieune et al.l (I1997LIT9 98~) colors (i.e., with¬ 
out color calibration). Compared to the calibrated isochrones, 
original model colors are generally bluer, and the difference 
becomes a few hundredths of a magnitude at the bottom of 
the MS (see bottom panels). In the VK S CMD, the best-fitting 
models tend to have bluer colors than the observed MS near 
MS turn-off (see top panels), which could be reconciled by an 
older age and/or a higher reddening. This is partly because our 
global \ 2 solution was heavily weighted by lower MS stars, 
in particular by stars in the Vic CMD with its longer color 
baseline than the other color indices. As noted above, we re¬ 
stricted fits to B - V < 1.1, V — Ic < 2.0, and V — K s < 2.2. 

5.6. Errors in the Cluster Parameters 

We took a 95.4% confidence interval in Figure[T9las a con¬ 
servative measure of a lrr fitting error in the cluster parame¬ 
ter estimates. These are shown in the first row (“fitting” er¬ 
ror) in Table [6] However, a total error budget for the clus¬ 
ter parameters is no longer dominated by fitting errors, by 
means of our empirical color corrections, but by the sum of 
various systematic errors such as those listed in the first col¬ 
umn of Table [6] They include systematic errors in the he¬ 
lium enrichment parameter (AT/AZ), reddening laws [By, 
Rvi = E(V - I C )/E(B - V), and R VK = E(V - K S )/E(B - V )], 
zero points in the cluster photometry [AV, A K s , A(B — V ), 
and A(V — Ic)], cluster parameters of the calibrating systems, 
and the bolometric correction of the Sun ( BCv,@ )• For alter¬ 
native values of each of these errors, we repeated the above 
search for the minimum \^ ot value using a downhill simplex 
method dPress et al 1119921) . and estimated a size of errors in 
[M/H], age, E(B - V), and ( m - M) o, as listed in the third 
through sixth columns. At the bottom of the table, we list 
total errors as a quadrature sum of these error contributions. 
Below we provide detailed information on each of these error 
sources. 

The helium abundance (T) sensitively affects isochrone lu¬ 
minosity by AMy /AT ~ 3 at fixed T e g- and heavy-element 
mass fraction (Z). To explore effects of helium abundance on 
our derived cluster parameters, we compared our MS fitting 
results with those using models with different helium enrich¬ 
ment parameters, AT/AZ = 2.0 and 3.5 (G. Newsham & D. 
M. Terndrup 2007, private communication). These models 
were constructed from the same YREC suite of programs, and 


many input parameters are similar to those used in our models, 
except that AT/AZ= 1.2 was assumed in our base models. At 
fixed T s s, the above models with higher AT / AZ are approxi¬ 
mately 0.1 and 0.3 mag fainter in bolometric magnitudes than 
our base models at [M/H]= +0.4. The a-element abundances 
in their models were set based on the abundance patterns ob¬ 
served for giants in the Galactic bulge dFulbright et al.1120071) . 
but [o/Fe] ratios at the metallicity of NGC 6791 are nearly 
solar. 

We applied the same empirical color-7)n corrections (Ta¬ 
ble 0 to the alternative models with AT/AZ = 2.0 and 3.5, 
and derived cluster parameters from the Xtot minimization. 
We obtained an error dependence by fitting a line to values 
of each cluster parameter as a function of AT /AZ. Errors 
in each quantity were then estimated assuming a(AY/AZ) = 
0.17, as derived from the primordial helium abundance, the 
initial solar value, and helium abundances of the Hyades and 
the Pleiades (see Paper III). A higher AT/AZ results in a 
fainter MS at a fixed color or a shorter distance to the clus¬ 
ter. However, the effect of helium on stellar isochrones is 
nearly monochromatic, meaning that photometric metallicity 
and reddening estimates are weakly dependent on AT/ AZ. 

The next three rows in Table [6] show effects of errors in the 
reddening laws. We adopted the same size of errors in Ry, 
Ryj , and Ryr as in Paper III. The next four rows list errors 
resulting from zero-point errors in the photometry. Within 
the fitting range (B — V < 1.10), we assumed 0.005 mag 
as a measure of systematic errors in the optical photome¬ 
try. ICarnev et al.l ( 120051) found that a weighted mean dif¬ 
ference between their CIT //-band photometry and 2MASS 
values transformed onto the CIT system (ICarpentedl200 lh is 
0.014 ±0.005 mag, with no apparent dependence on magni¬ 
tude. This difference is larger than the calibration uncertainty 
(0.007 mag) that was specified in the explanatory supplement 
to the 2MASS All Sky Data Release[|j We took an error in 
K s as a quadrature sum of these two errors. The sense in Ta¬ 
ble [6] is that if B — V photometry were systematically under¬ 
estimated by 0.005 mag (bluer) than true values, a true metal¬ 
licity is 0.028 dex higher than the photometric metallicity. 

Calibration errors in Table [6] are a quadrature sum of pa¬ 
rameter errors from calibrating cluster systems employed in 
this study. Since models matched to NGC 6791 CMDs were 
primarily calibrated using both the Hyades and Praesepe, we 
took errors from these two clusters: for the Hyades, we 
adopted errors of <r([Fe/H]) = 0.01, a[E(B — F)] = 0.002, 
a[(m — M) o] = 0.01. For Praesepe, we adopted <r([Fe/H]) = 
0.02, a[E(B - V)] = 0.002, a[(m - M) 0 ] = 0.04. We simply 
took 100 Myr for an error in age from the consideration that 
both clusters are about 550-650 Myr old. 

There can be a zero-point shift in the absolute V magni¬ 
tude scale, depending on the adopted absolute visual magni¬ 
tude of the Sun. The My of the Sun in our solar-metallicity 

10 See http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html 
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Fig. 20.— Isochrone fits on CMDs of NGC 6791. Only those left from the photometric filtering routine are shown with error bars, representing a single¬ 
star cluster sequence. Solid red lines are calibrated isochrones with the best-fitting parameters (without binary corrections): [M/H]= +0.38, E(B — V) = 0.108, 
(m — M )o = 13.00, and the age of 9.5 Gyr. The sam e models with ±0.5 Gyr in age are additionally shown by solid red lines in the top panels. The dotted-dashed 
blue lines are the sam e theoretical models, but with ILeieune et al] <19971Il998l) colors. The diagonal dotted lines on the right panels represent the completeness 
limit in K in ICamey et al.l <20051) . 


TABLE 6 

Error Budget for the Cluster Parameters of NGC 6791 


Source of 
Error 

Adopted 
Size of Error 

A[M/H] 

(dex) 

AAge 

(Gyr) 

A E(B-V) 
(mag) 

A(m—M)o 

(mag) 

Fitting.... 


±0.015 

±0.20 

±0.004 

±0.023 

AY/AZ.. 

±0.17 

±0.001 

±0.02 

±0.002 

4=0.016 

Rv . 

±0.30 

±0.022 

4=0.07 

4=0.007 

4=0.006 

Rvi . 

±0.07 

±0.027 

4=0.16 

4=0.007 

±0.025 

Rvk . 

±0.12 

±0.002 

±0.07 

4=0.002 

4=0.002 

AV . 

±0.005 

4=0.002 

4=0.02 

±0.001 

±0.005 

A (B-V). 

±0.005 

±0.028 

4=0.04 

4=0.004 

±0.004 

A (V-I)c- 

±0.005 

4=0.021 

±0.10 

±0.006 

4=0.018 

A K s . 

±0.016 

±0.005 

±0.06 

4=0.002 

±0.001 

Calibration 


±0.022 

±0.10 

±0.003 

±0.041 

BCV,0 .... 

±0.04 




±0.040 

Total . 


±0.056 

±0.32 

±0.014 

±0.071 
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model is My t 0 = 4.86 at = 5777 K assuming a solar 
age of 4.57 Gyr. This value is 0.05 mag larger than My,o = 
4.81 ±0.03 in iTorres! (120101) . We adopted the difference as 
a measure of a systematic error in the bolometric correction 
(BCy ,q). and included it in Table [6] 

In addition to the above systematic errors, unresolved clus¬ 
ter binaries and/or photometric blends can modify MS-fitting 
parameter estimates, as they are brighter than single MS stars 
at a given color, forming a skewed distribution above a sin¬ 
gle MS. Although our photometric filtering routine reduced a 
number of blended sources on CMDs, low mass-ratio bina¬ 
ries in particular could remain and affect our parameter esti¬ 
mation. iBedin et aD (i2008l) estimated a binary fractiorjT] of 
NGC 6791 along the MS, and found 32 ±3% in the core 
of the cluster. Because the fraction is rather uncertain, we 
adopted a 40% binary fraction like those used for nearby sys¬ 
tems in Paper III. At a 40% binary fraction, results from arti¬ 
ficial cluster tests (Paper III) showed that photometric metal- 
licities are underestimated by ~ 0.036 dex, while E(B - V) 
are overestimated by ~ 0.003 mag even after the photometric 
filtering. The corresponding change in distance modulus is 
A(m - AT) o « +0.04 mag. The effect on age was not included 
in the simulation, but is likely to be small. 

Together with the calibration and the bolometric correction 
errors in Table [6] the binary correction is one of the major 
systematic uncertainties in the determination of the cluster’s 
distance. However, we caution that binary corrections are 
uncertain owing to the poorly constrained nature of binaries 
in the cluster, such as a binary fraction, a mass function for 
secondaries, and a distribution of binary mass ratios. Nev¬ 
ertheless, the above biases are systematic in nature, and we 
list adjusted values of the cluster’s metallicity, reddening, and 
distance with binary corrections in the second row of Table [5] 
We added in quadrature the size of the binary corrections to 
the total error budget. The cluster CMDs are relatively free of 
foreground/background stars unlike the visual impression of 
CMDs. We contend that effects of these non-cluster members 
are small compared to the other systematic errors. 

6. SUMMARY AND DISCUSSIONS 

We extended our effort to calibrate stellar isochrones in the 
Johnson-Cousins and 2MASS filter systems, using a set of 
photometric data for well-studied open cluster systems. Our 
base models were constructed using the YREC suite of pro¬ 
grams, where model T c tj were initially converted into observ- 
able colors using color-T e ff transformations in iLeieune et al.l 
(1997. 1998 ). In our earlier series of papers, we defined color 
corrections as a function of 7/if to match photometry of the 
Hyades and the Pleiades. In this work, we added Praesepe, 
extending earlier color calibration to cooler MS stars. We 
focused our analysis on the test of the models, and thus our 
empirical color corrections, using cool and metal-rich stars in 
NGC 6791. 

There exists no geometric distance measurement for 
NGC 6791. As an alternative, we utilized an accurate metal¬ 
licity of the cluster from a number of recent high-resolution 
spectroscopic studies to validate our calibrated models. Since 
optical and near-IR color indices have different sensitivi¬ 
ties to metallicity, we examined whether the cluster’s dis¬ 
tance and reddening estimates from different CMDs converge 
on each other at the spectroscopic metallicity of the cluster. 

11 The binary fraction is defined as the number of binaries divided by the 
total number of systems. 


When the original ILeieune et al.l (119971119981) color-T e ff rela¬ 
tions were used to transform theoretical T e ff into colors, we 
found internally inconsistent results in distance and redden¬ 
ing. However, we demonstrated that our empirical color cor¬ 
rections to the ILeieune et all relations bring about the con¬ 
sistency, improving the accuracy of MS-fitting distance and 
reddening estimates. We also found that YREC models are 
in satisfactory agreement with the observed mass-radius and 
mass-luminosity relations from the two binary systems in 
NGC 6791, when reasonable assumptions on T e ff scale and 
helium abundance were made. The above tests support our 
earlier conclusion that discrepancies found between models 
and observed CMDs mainly originate from systematic errors 
in the color-r e ff transformations. Our calibrated isochrones 
are valid at least within a metallicity range covered by the cal¬ 
ibrating cluster systems and NGC 6791 (0 < [Fe/H] < +0.4). 

The accuracy of MS fitting can be improved by our empir¬ 
ical color-T e ff corrections. In particular, differential sensitiv¬ 
ities of broadband colors enabled us to constrain reddening, 
distance, and metallicity of NGC 6791. Our photometric so¬ 
lution yields E(B - V) = 0.105 ± 0.014, [M/H]= +0.42 ± 0.07, 
( m - M) o = 13.04 ± 0.09, and the age of 9.5 ± 0.3 Gyr, after 
excluding B — V colors of the cool MS stars in the cluster. 
Errors in these quantities include all potential systematic er¬ 
rors involved in the MS fitting. Our best-fitting age may be 
~ 1 Gyr too old since our models do not incorporate micro¬ 
scopic diffusion. Our E(B — V ) estimate is found at the lower 
boundary of earlier E(B — V) estimates in the literature, but 
the lower reddening also leads to a better agreement of mod¬ 
els with the eclipsing binary data, when temperature s are es - 
timated from the IRFM relations in lCasagrande et all d2010t) . 

However, we found that the best-fitting models from the up¬ 
per MS do not match well observed B — V colors in the lower 
MS of NGC 6791. We could not fully resolve this unfortu¬ 
nate circumstance. Nevertheless, an excellent agreement in 
V — Ic for the same set of stars, as well as tests using IRFM 
color-T e ff relations and independent cluster photometry, sug¬ 
gests that the problem could originate from a color-scale error 
in the lStetson et al.l photometry. While the photometric differ¬ 
ence of ~ 0.02 mag could still be considered to be acceptable 
in many applications, our empirical corrections are probably 
better defined. For example, the B - V correction in the top 
panel of Figure QT| shows that the scatter of individual points 
around the mean line is ~ 0.016 mag at 7)n ~ 4500 K. If our 
adopted cluster photometry is correct, the dispersion trans¬ 
lates into an error in the mean of 0.007 mag in a 4-8 point 
moving window. 

If stars with ~ 0.02 mag photometric offsets in B — V were 
included in MS fitting, our photometric solution of the clus¬ 
ter parameters would yield a smaller metallicity (A[M/H] ~ 
0.1 dex), a larger foreground reddening (A E(B — V) ~ 0.01), 
and a shorter distance modulus (A (m - M) o ~ 0.1), accord¬ 
ing to our error matrix in Table [6] Since the presence of such 
photometric offsets could reduce the accuracy of MS fitting, 
we excluded B - V colors of the red MS stars in the above 
parameter estimates. Until more evidence is collected, our 
calibrated isochrones can be used with limited accuracy for 
the cool and metal-rich MS stars. 

Our best-fitting distance modulus for NGC 6791 is in good 
agreement with lBrogaard et ak i (1 2011 ). who found (m - M) o = 
13.01 ±0.08 from eclipsing binary systems. They inferred 
luminosities of individual eclipsing binary components from 
measured stellar radii and spectroscopic temperatures, and 
adopted bolometric corrections in V from MARCS models to 
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derive absolute V magnitudes of stars. In their derivation of 
the cluster’s distance, they assumed E(B - V) = 0.160 ±0.025 
from a comparison between spectroscopic and photometric 
temperatur es of the binary compo n ents. Based on the dista nce 
estimate in Brogaard et al.l (1201 11) . IBrogaard et al.1 (120121) ex¬ 
tended their analysis of NGC 6791 by combining the clus¬ 
ter’s CMDs. They employed Victoria-Regina and DSEP 
isochrones jDotter et al.1120081) . and found that the cluster 
is about 8.3 Gyr old, and has a helium mass fraction of 
Y = 0.30 ±0.01 when the metallicity is [Fe/H]= +0.35 with 
[Grevesse & Sauval (119981) abundance mixtures. From the 
isochrone fitting on the Vic CMD, they derived E(B -V) = 
0.14 ± 0.02 for the global mean of the color excess^_ 

More recent work (Casagran de & VandenBer d 1201 4t 
IVandenBerg et al.1 1201 4l) claimed that equally acceptable 
isochrone fits can be obtained with slightly different pa¬ 
rameter inputs from those in iBrogaard et ali (1 2012 1. al¬ 
though they did not fully explore all _available parame¬ 
ter space. For example, ICasagrande & VandenBergl d2014l) 
found that an 8.5 Gyr old Victoria-Regina_model at Y = 
0.28 and [Fe/H]= +0.30 with lAsplund et aU (12009b abundance 
mixtures shows nearly indistinguishable fits when E(B — 
V) = 0.16 is assumed o n BVI J and ugriz CMDs. Simi¬ 
larly, VandenBerg et al. (2 014) used the same model as in 
Casag rande & V andenBerg to find a good fit with E(B -V) = 
0.14 from the VJ CMD, and attr ibuted the difference in 
E(B-V) from lCasagrande & VandenBergl to a potential zero- 
point error in /-band photometry. 

These recent E(B — V) estimates, although not statistically 
significant given the size of our estimated total errors, are sys¬ 
tematically higher than our best-fitting solution [E{B - V) = 
0.105 ±0.0141 and other previous measurements (IJanesll984t 
Montgomery et al.1 1994t IChabover et al.Jll999t IStetson et al.1 
20031:ICarraro et al 1120061: iMontalto et al.fl2007l) . The ~ 0.03- 
0.05 mag difference in E(B - V ) is partly due to different 
model isochrone s employed in each of th e studies. Although 
Victoria-Regina (IVandenBerg et al.1120141) and YREC mod¬ 
els agree with each other on the bolometric magnitude ver¬ 
sus temperature plane at the metallicity of NGC 6791 (as¬ 
suming the same helium abundance), the former models with 
color-/)n conversions from MARCS produce slightly bluer 
(~ 0.02 mag) and brighter magnitudes (~ 0.1 mag) near the 
cluster’s MS turn-off in the BVl ( CMDs. 

However, it should be emphasized that our empirical color 
corrections were defined to bring the isochrones to match 
the observed MS of nearby clusters, instead of relying solely 


on theoretical computations (e.g., Casagrande & VandenBergl 
120141) . and therefore should provide more reliable pre¬ 
dictions on colors and magnitudes of stars. In addi- 
tion, we note that isochrone fitt i ng in the above studies 
(ICasagrande & VandenBergl 12014 : IVandenBerg et al.1 120141) 
focused on relatively bright stars in the cluster (V < 20 or 
T e ff > 4700 K in Vic CMD), while their best-fitting isochrones 
show systematic departures for cooler stars. On the other 
hand, our work in this paper concerns model matches down 
to V - Ic ~ 2.0 or V ~ 22 mag. 

Distances estimated using our calibrated models are ulti¬ 
mately tied to the Hipparcos distance to the Hyades and the 
mean geometric distance to the Pleiades (see Paper III). The 
distance to Praesepe adopted in the current work was derived 
using the Hyades-based calibration. Our luminosity calibra¬ 
tion helps not only to determine distances to clusters and in¬ 
dividual stars, but also to make an accurate estimate of their 
foreground reddening. Both quantities are necessary to have 
a proper luminosity calibration of astrophysical objects. Nev¬ 
ertheless, there is an important limitation in our methodol¬ 
ogy that the current empirical corrections rely on a set of 
calibrating clusters near solar metallicity ([Fe/H]« +0.04 or 
+0.14). Although we demonstrated the accuracy of our mod¬ 
els at [Fe/H]~ +0.4 using different color indices, our calibra¬ 
tions should ultimately be tested and confirmed using metal- 
rich stars with accurately known distance and reddening. Cur¬ 
rently, there are only a few star clusters that can be used for ac¬ 
curate calibrations of stellar colors and magnitudes with pre¬ 
cise cluster parameters. However, there will be more cluster 
systems available in the era of Gaia dPerrvman et al.1[2001 ). 
over a wider range of stellar mass, metallicity, and age. 
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